{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "205909dd",
   "metadata": {},
   "source": [
    "#  第五讲 无标度网络【实践】二"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d43cada7",
   "metadata": {},
   "outputs": [],
   "source": [
    "#!/usr/bin/python3\n",
    "# -*- coding: utf-8 -*-\n",
    "# Author ： 单哥的科研日常\n",
    "# 关注B站和公众号：单哥的科研日常，获取更多讲解教程"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "6f9db701",
   "metadata": {},
   "source": [
    "### 实验环境：\n",
    "### Python版本==3.9.6, networkx==2.6.3, matplotlib==3.5.2, numpy==1.23.1, pandas==1.4.3, scipy==1.9.0"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "baadd879",
   "metadata": {},
   "source": [
    "## 1、配置模型生成给定度序列的网络"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "1752884b",
   "metadata": {},
   "outputs": [],
   "source": [
    "# 导入库\n",
    "import networkx as nx\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "from networkx.utils import powerlaw_sequence\n",
    "from utils import *"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "61fe2885",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(641, 641)\n",
      "(1006, 1006)\n",
      "(1006, 1006)\n",
      "(1006, 1006)\n",
      "(1006, 1006)\n",
      "(1138, 1138)\n",
      "(1511, 1511)\n",
      "(1725, 1725)\n",
      "(1725, 1725)\n",
      "(1725, 1725)\n",
      "(1725, 1725)\n",
      "(1725, 1725)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1780, 1780)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(1824, 1824)\n",
      "(2312, 2312)\n",
      "(2996, 2996)\n",
      "(3173, 3173)\n",
      "(3479, 3479)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4010, 4010)\n",
      "(4965, 4965)\n"
     ]
    }
   ],
   "source": [
    "n, gamma = 5000, 2.1\n",
    "degree_seq = powerlaw_sequence(n, gamma)\n",
    "int_deg = [int(di) for di in degree_seq]\n",
    "\n",
    "while sum(int_deg)%2 != 0:\n",
    "    degree_seq = powerlaw_sequence(n, gamma)\n",
    "    int_deg = [int(di) for di in degree_seq]\n",
    "    \n",
    "    \n",
    "G = nx.configuration_model(int_deg)\n",
    "# 该模型会生成自环和多重链接，而这些在真实网络中通常是不会出现的\n",
    "for e in G.edges():\n",
    "    if e[0]==e[1]:\n",
    "        print(e)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "f0665bbb",
   "metadata": {},
   "outputs": [],
   "source": [
    "k, Pk = get_pdf(G)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "1c6a6f33",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAE4CAYAAABIYsQpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAXN0lEQVR4nO3db2gl133G8eenu5KJbBCNZPrCtiSHNSEm4KZWE9pCyYv8cdxsEtKS2qiBehcrXpPgvqtTBdpCt25oXyQmbhwl3jrpqg7GlGbdbjC0YExat6ycvKhT4+xiLFsYss66iCZLsXf31xeja11dzdw59+7c+XPm+4HL3XvuuTNnjhc/O3PmnDF3FwAAeSaqbgAAoBkIDABAEAIDABCEwAAABCEwAABBCAwAQJADVTdgnObm5nxxcbHqZgBAozz33HM/c/dr+8ujDozFxUVtbGxU3QwAaBQz20wr55IUACAIgQEACBJlYJjZITNb297erropABCNKAPD3Z9095WZmZmqmwIA0YgyMAAAxSMwAABBCAwAQJAoA+OKBr3X16XFRWliInlfXy+6eQDQSFEGxsiD3uvr0sqKtLkpuSfvKyuEBgAo0sAY2eqqdOHC3rILF5JyAGg5AqPXK68MVw4ALUJg9JqfH64cAFokysAYedD72DFpenpv2fR0Ug4ALRdlYIw86L28LK2tSQsLklnyvraWlANAy0W9vPlIlpcJCABIEeUZBgCgeAQGACAIgQEACEJgAACCRBkYPEAJAIoXZWDwACUAKF6UgQEAKB6BAQAIQmAAAIIQGACAIAQGACAIgQEACBJlYDAPAwCKF2VgMA8DAIoXZWAAAIpHYAAAghAYAIAgBAYAIAiBAQAIQmAAAIIQGACAIAQGACAIgQEACEJgAACCEBgAgCBRBgaLDwJA8aIMDBYfBIDiRRkYtbe+Li0uShMTyfv6etUtAoBcB6puQOusr0srK9KFC8nnzc3ksyQtL1fXLgDIwRlG2VZXd8Oi68KFpBwAaozAKNsrrwxXDgA1QWCUbX5+uHIAqAkCo2zHjknT03vLpqeTcgCoMQKjbMvL0tqatLAgmSXva2sMeAOoPe6SqsLyMgEBoHE4wwAABCEwAABBCAwAQBACAwAQhMAAAAQhMAAAQRoTGGb2LjN7xMyeqLotANBGpQSGmR03s3Nm9nxf+W1m9qKZnTWz+wdtw91fcvcj420pACBLWRP3HpX0NUnf6RaYWUfSQ5I+LGlL0mkzOympI+mBvt8fdvdz5TQVAJCmlMBw92fMbLGv+P2Szrr7S5JkZt+V9El3f0DSx8toFwAgXJVjGNdJerXn89ZOWSozmzWzhyW9z8y+OKDeipltmNnG66+/XlxrAaDlqlxLylLKPKuyu5+XdE/eRt19TdKaJC0tLWVuDwAwnCrPMLYk3dDz+XpJr1XUFgBAjioD47Skm8zsRjObknSHpJMVtgcAMEBZt9U+JulZSe82sy0zO+LuFyV9XtJTkl6Q9Li7/7ig/R0ys7Xt7e0iNlcf6+vS4qI0MZG8r69X3SIALWLu8V7mX1pa8o2NjaqbUYz1dWllRbpwYbdsepqHLwEonJk95+5L/eWNmendequre8NCSj6vrlbTHgCtE2VgRHlJ6pVXhisHgIJFGRju/qS7r8zMzFTdlOLMzw9XDgAFizIwonTsWDJm0Wt6OikHgBIQGE2xvJwMcC8sSGbJOwPeAEoU5V1SZnZI0qGDBw/efebMmaqbAwCN0qq7pKIcwwCAikUZGACA4hEYAIAgBAYAIEiUgRHlxD0AqFiUgcGgNwAUL8rAAAAUj8CIFUuhAyhYlY9oxbj0L4W+uZl8lpgZDmBkUZ5htH7Qm6XQAYxBlIHR+kFvlkIHMAZRBkbrsRQ6gDEgMGLEUugAxoDAiBFLoQMYA+6SitXyMgEBoFCcYQAAgkQZGK2/rRYAxiDKwGj9bbUAMAZRBgYAoHgERlux1hSAIXGXVBux1hSAEXCG0UasNQVgBARGG7HWFIAREBhtxFpTAEZAYLQRa00BGEGUgcHEvRysNQVgBObuVbdhbJaWlnxjY6PqZgBAo5jZc+6+1F8e5RkGAKB4BAYAIAiBAQAIQmAAAIIQGACAIAQGACAIgQEACEJgYC+WPQeQgeXNsYtlzwEMwBkGdrHsOYABogwM1pIaEcueAxggysBw9yfdfWVmZqbqpjQLy54DGCDKwMCIWPYcwAAEBnax7DmAAbhLCnstLxMQAFJxhgEACEJgIAwT+oDW45IU8jGhD4A4w0AIJvQBEIGBEEzoAyACAyGY0AdABAZCpE3om5yUfv5zBsGBFiEwkK9/Qt/sbPJ+/rzkvjsITmgAUSMwEGZ5WXr5ZenyZemaa6Q339z7PYPgQPQIDAyPQXCglYYKDDObGldD0CAMggOtNOwZxp91/2Bmv1lwW9AUrGoLtNKwgfGUma2Y2bKkj42jQVnM7FNm9k0z+56ZfaTMfaMPq9oCrWTuHlbR7BFJ25J+RdJ/uPsfB+/E7Likj0s65+7v7Sm/TdJXJXUkfcvd/zJgW78k6a/d/Uhe3aWlJd/Y2AhtJgBAkpk95+5L/eXBa0m5+xEze4ekX5X0a2b2DXf/XODPH5X0NUnf6WlQR9JDkj4saUvSaTM7qSQ8Huj7/WF3P7fz5y/t/A4AUKLcS1Jm9jEz+08ze1HStyVddvevDBEWcvdnJL3RV/x+SWfd/SV3f1PSdyV90t3/y90/3vc6Z4kvS/q+u/8w/BBRGla0BaIWcobxN5J+X9J/S7pV0l+Z2UPu/tgV7vs6Sa/2fN6S9IEB9b8g6UOSZszsoLs/nFbJzFYkrUjSPHftlIcVbYHohQx6/9Td/83d/8fd/0XSRyUVMUPLUsoyB1Tc/UF3v9Xd78kKi516a+6+5O5L1157bQHNRBBWtAWiFxIYL5vZn/fMwXhL0v8WsO8tSTf0fL5e0msFbBdVYDIfEL2QwHBJn5b0qpn9QNJZSU+b2U1XuO/Tkm4ysxt3wugOSSevcJuSJDM7ZGZr29vbRWwOIZjMB0QvNzDc/U53v1nSgqQ/VDJ572pJ3zKzVwf9tsvMHpP0rKR3m9mWmR1x94uSPi/pKUkvSHrc3X882mHsa/OT7r4yMzNTxOYQIm0yn1kylsEAOBCFYW6r/T9JGzuvobj7nRnlpySdGnZ7qKHuwPbqahISZslKthID4EAkolx8kEtSFemuaLuwsBsWXQyAA40XZWBwSapiDIADUYoyMFAxBsCBKBEYKB6r2QJRijIwGMOoGKvZAlEKXq22iVitFgCGl7VabZRnGACA4hEYKAcr2QKNFzxxDxgZK9kCUYjyDINB75phJVsgClEGBhP3aoaJfEAUogwM1MygiXzdsQ0z6cCB5J0xDqCWCAyMX9ZEvttvT8YyNjeTskuXkvfuGAehAdQKgYHxy5rId+rU/rGNLsY4gNqJcuKemR2SdOjgwYN3nzlzpurmIMvExP5VbXuZSZcvl9ceAJJaNnGPQe+GyFuMkMUKgVqJMjDQEGljG10sVgjUDoGB6vSObUhSp5O8s1ghUEsEBsrXu0zI6mpyJuEuXbyYvB87lpSzjAhQKywNgnLlLRPCMiJAbUV5l1QXy5vX0OLi7ryLXgsLyfPA874HMHatukuKtaRqLG+ZkKzvNze5NAVULMrA4LbaGst73vegW2mZ/Q1UKsrAQI3lPe970K22zP4GKkVgoFx5z/vufp+FFW6ByjDojXpi8BuoTKsGvRGBvEtXAEpHYKCe8i5dASgdE/dQX8vLBARQI5xhAACCRBkYTNwDgOJFGRhM3AOA4kUZGACA4hEYaJ7e5dGzlj8PqQNgKNwlhWYJWf6cJdKBsWCmN5olZAY4s8SBK8JMbzRX7+WltCCQ9q4xlbeEOoCREBiot+7lpc3N5PGtWXqXRc9bQh3ASAgM1Nvq6u5YRJb+NaZYhwoYCwID9TboMlLWGlOsQwWMBXdJod7m50cbwGYdKqBwnGGg3ri8BNRGlIHBWlIRCbm81HsX1dxc8mLCHlA45mGg2fon6fWbnmb8AhgS8zAQp7y7qC5cSOoAuGIEBpotZDIeE/aAQhAYaLaQyXiD6vSPf1xzTTJWYpZ8ZgwEeBuBgWZLu4uq16A7qvpnkZ8/L/3iF7vfnz8vHT5MaAA7CAw0W/9dVLOzyStkwl7ILPI332QMBNjBXVJor4mJwetTdZlJly+Pvz1ATXCXFOKTN/9ifT0pyxqTCF2MkEULAUksDYKm6p9/cf787nebm9JddyVnBZcu7ZZ3xySk5DLVsWOD53BI0tQUs8qBHZxhoJnyxh/eemtvWHT1jkmkjX9cffVu3dlZ6fhxJv0BOxjDQDOFjj+kYUwCGIgxDDRf75jFxBX81e3+Pmutqd79sB4V8DbGMNAM/WMWaZebek1O7h/D6OqWbW4m25R2Lzv17yetDtBSXJJCMywupj8Xo9NJguGd70w+v/FGcldTd6D6vvt2B8QnJtIvRfU+WyNrP3nP3wAiknVJijMMNEPWelCXLw8ej+g9K8i6jNW77az9sB4VwBgGGiJrLsQwcyRCtlHEfoBINSYwzOw9ZvawmT1hZkerbg9KVsST90K2wRP+gEylBIaZHTezc2b2fF/5bWb2opmdNbP7B23D3V9w93skfUbSvmtriFzIk/eG3cbsrPSOd0if/ezu3VBF7AeIVCmD3mb2W5J+Luk77v7enbKOpJ9I+rCkLUmnJd0pqSPpgb5NHHb3c2b2CUn3S/qau/993n4Z9EamtCf18XQ+QFL2oHdpd0mZ2aKkf+oJjF+X9Kfu/tGdz1+UJHfvD4u0bf2zu/92Xj0CA5m4GwrIVMe7pK6T9GrP5y1JH8iqbGYflPRpSVdJOjWg3oqkFUmaZ6ASWbgbChhalYFhKWWZpzvu/rSkp/M26u5rktak5AxjxLYhdvPz6WcY/CMDyFTlXVJbkm7o+Xy9pNcqagvahruhgKFVGRinJd1kZjea2ZSkOySdLGLDZnbIzNa2t7eL2BxiNOzdUIPWl2LtKbREWXdJPSbpg5LmJP1U0p+4+yNmdrukryi5M+q4uxf6zzsGvVGIQXdUSdxthehUfpdUFQgMFGLQHVUSd1shOnW8S2pszOyQpEMHDx6suimIwSh3VHG3FSLUmKVBhuHuT7r7yszMTNVNQQwGrS/F2lNokSgDAyjUoDuquNsKLRLlJSmgUN3B69XV5FJT93kbvYPag74DIhHloHfPGMbdZ86cqbo5ANAorXqmN2MYAFC8KAMDSDXM5Lt77x1cd24umfDX+5qbY9IeosYYBtqhf/Ld5mbyuav/u69/ffe7/rp33SW99db+fZw/Lx0+nPyZMQxEKMoxjC4m7uFto0y+G7Uuk/bQcEzcQ7sVsZx5aF0m7SFSUY5hMOiNfUaZfDdqXSbtIVJRBgawz7CT7/r11p2czK43NcWkPUQryktSwD7DTr67/Xbp1Knsuvfdlwxy95qdlb76VQa8ES0GvQEAe7Rq4h4AoHhRBgZP3AOA4kUZGNwlBQDFizIwAADFIzCAQWtM5dXpLZ+bk666andtqU4nWZOqqH2bJfXqtnZVyDEgDu4e7evWW291YKATJ9ynp92l3df0dFKeV+fo0f3laa+jR4vdd+9rampv/bKFHAMaR9KGp/w/ldtq0W6D1pjqrgeVVafTkS5dyt9HpyNdvFjsvrPqly3kGNA4WbfVRhkYPEAJwSYmkn8X9zOTLl8eXGcYab8vat+99csWcgxonFbNw3DukkKoQWtM5dXpdML2kVXvSvYdsp0yhBwDohFlYADBBq0xlVdnZSV/DSpp77M0ith3r6rXrgo5BsQjbWAjlheD3ghy4oT7woK7WfKeNmCbVae3fHY2GYTuDv5OTGQPeI+ybymp193+7Gw9BpdDjgGNIga9AQAhWjWGAQAoHoEBAAhCYAAAghAYAIAgUQYGy5sDQPGiDAxn4h4AFC7KwAAAFI/AAAAEITAAAEEIDABAEAIDABCEwAAABCEwAABBCAy0w/p68jjRiYnk/d57935eX8/+3dxc8gQ5s+TPg+ouLib1JiZ2f9N9dTr7y9Je3XoHDuzus9uGblnaq3scacfaewzXXJO80o4p73h7t91tV/fPvdvM6qv+tmX1Zdl6/9t1+7hO7evK679x92/amuexvHgeBtw9eT7D9PTucyTSXtPT+5/jcOKE++Tk/rpTU+l18/ZRxmtycu8zOUJfU1PJszsGHe8ox9jbV2m/T+v3sg06rjq0ryuv/wrsX/E8DLTW4qK0uZlfb2FBevnlsN8NU7cpOh3p0qX07xYWkvdRjrHbV1l91N+XZcv7b1d1+7ry+q/A/s16HkaUgWFmhyQdOnjw4N1nzpypujmo2sRE8u+tPGbS5cthvxumbgzMkvdRjrHbV1l91N+XZcv7b1d1+7ry+q/A/m3VA5SctaTQa35+tHqDfjdM3abodLK/m58f/Ri7v8v6fdV9l7f/qtvXldd/JfRvlIEB7HHsmDQ9PbjO9HRSr/93k5P7605NpdfN20cZJieT9g1rakpaWRl8vKMcY29fpf0+rd/LNui46tC+rrz+K6N/0wY2Ynkx6I23nTjhvrDgbpa8Hz2693PWwOCJE+6zs7uDiLOzg+suLCT1zPYPoE5MhA0Ud+t1Orv77LahW5b26h5H2rH2HsPVVyevtGPKO97ebXfb1f1z7zaz+qq/bXUaUO7+t+v2cZ3a15XXfwX1rxj0BgCEaNUYBgCgeAQGACAIgQEACEJgAACCEBgAgCAEBgAgCIEBAAgS9TwMM3td0qakGUnbfV/3l/V/npP0s7E2ML1dRf8ur27W96HleX1b534c5rfj7se0srL/Tta5H7O+q2M/ZrWryN+Nux8X3P3afb9Mm80X20vSWl5ZyufUmY7jblfRv8urm/V9aHle39a5H4f57bj7MaRvx92Xde7H0D6rQz9eSV/WqR/TXm25JPVkQFlanXEbdZ/D/C6vbtb3oeUhfTtuV7K/0N+Oux/TyujH/O/q2I9Xss869eM+UV+SuhJmtuEpU+MxHPqxOPRlMejH0bXlDGMUa1U3IBL0Y3Hoy2LQjyPiDAMAEIQzDABAEAIDABCEwAAABCEwApnZ1Wb2bTP7ppktV92epjKzd5nZI2b2RNVtaTIz+9TO38XvmdlHqm5PU5nZe8zsYTN7wsyOVt2eumt1YJjZcTM7Z2bP95XfZmYvmtlZM7t/p/jTkp5w97slfaL0xtbYMP3o7i+5+5FqWlpvQ/bjP+78XfwDSb9XQXNra8h+fMHd75H0GUncapuj1YEh6VFJt/UWmFlH0kOSPibpZkl3mtnNkq6X9OpOtUsltrEJHlV4PyLboxq+H7+08z12Paoh+tHMPiHpB5L+tdxmNk+rA8Pdn5H0Rl/x+yWd3fmX8JuSvivpk5K2lISG1PJ+6zdkPyLDMP1oiS9L+r67/7DsttbZsH8f3f2ku/+GJC415+B/fPtdp90zCSkJiusk/YOk3zGzr6uapQaaJrUfzWzWzB6W9D4z+2I1TWuUrL+PX5D0IUm/a2b3VNGwhsn6+/hBM3vQzL4h6VQ1TWuOA1U3oIYspczd/ReS7iq7MQ2W1Y/nJfE/uHBZ/figpAfLbkyDZfXj05KeLrcpzcUZxn5bkm7o+Xy9pNcqakuT0Y/FoB+LQT8WgMDY77Skm8zsRjObknSHpJMVt6mJ6Mdi0I/FoB8L0OrAMLPHJD0r6d1mtmVmR9z9oqTPS3pK0guSHnf3H1fZzrqjH4tBPxaDfhwfFh8EAARp9RkGACAcgQEACEJgAACCEBgAgCAEBgAgCIEBAAhCYAAAghAYAIAgBAZQIjP7kJn9XdXtAEZBYADlukXSj6puBDAKAgMo1y2SfmRmV5nZo2b2F2aWtvQ2UDs8DwMo1y2SzilZBO9b7n6i4vYAwVh8ECiJmU1K+pmkTUmfc/dnK24SMBQuSQHluVnJcxkuSrpUcVuAoREYQHlukfTvSh7e87dm9ssVtwcYCoEBlOcWSc+7+08k/ZGkx3cuUwGNwBgGACAIZxgAgCAEBgAgCIEBAAhCYAAAghAYAIAgBAYAIAiBAQAIQmAAAIL8P9iNvhbp+v8lAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x345.6 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(6,4.8))\n",
    "plt.plot(k, Pk, 'ro')\n",
    "plt.xlabel(\"$k$\")\n",
    "plt.ylabel(\"$p_k$\")\n",
    "plt.xscale(\"log\")\n",
    "plt.yscale(\"log\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e664895a-e3c4-4cfe-ac36-174a457fff23",
   "metadata": {},
   "source": [
    "### **总结**：配置模型会生成自环和重连边，而这些在真实网络中通常是不会出现的。生成的网络的平均度会随着幂律指数的改变而变化，即平均度不可控"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "03c40215-f1c1-4150-b687-71ea453aa908",
   "metadata": {},
   "source": [
    "## 2、隐参数模型生成给定度分布指数的无标度网络"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "fd667cb4-851e-4bf5-9d0d-a57e7f7d5130",
   "metadata": {},
   "outputs": [],
   "source": [
    "N = 5000\n",
    "gamma = 2.1\n",
    "avk = 6.0\n",
    "L = int(avk*N/2)\n",
    "\n",
    "G = generate_SF_network(N, gamma, L)\n",
    "k, Pk = get_pdf(G)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "7e542219-ea2f-4183-be20-27a8a4fd80fa",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAE4CAYAAABIYsQpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAW9klEQVR4nO3db4wc9X3H8c/31ndpDqR7cJf2AeBbkBEKakUjTqhqHzQPkhSiGKIUJUQrRILFxUakD/qkpE4fVKqLoj6pUHDIpdgmZgtCqEpMS4SUSAilpRXn5EFBiNpCPnNCqo1TnUhJZfB9+2BuuL29mb3f3M7OzM68X9Jqvb+d/c3PN5w/zP7+mbsLAICdTJTdAADAeCAwAABBCAwAQBACAwAQhMAAAAQhMAAAQfaU3YBRmpub83a7XXYzAGCsnD59+l13/0R/ea0Do91ua3l5uexmAMBYMbOVpPJafiVlZvvNbGltba3spgBAbdQyMNz9eXdfnJmZKbspAFAbtQwMAED+CAwAQBACAwAQhMAAAAQhMAAAQQiMft2u1G5LExPRc7dbdosAoBIIjF7drrS4KK2sSO7R8+JiMaFBUAGoOAKj1+HD0vvvby17//2ofJTKDCoACERg9Dp/Plt5XsoKKgDIgMDotXdvtvK8lBVUAJABgdHryBFpenpr2fR0VD5KZQUVAGRAYPTqdKSlJWl+XjKLnpeWovJRKiuoACADc/ey25A7M9svaf++ffseOHPmTNnNCdPtRn0W589HdxZHjow+qAAggZmddveFbeV1DIzYwsKCsx8GAGSTFhh8JQUACEJg1AmT/wCMUK23aG2UePJfPJ8jnvwn0RcCIBfcYdQFk/8AjBiBURdM/gMwYgRGXTD5D8CIERh1weQ/ACNGYNRFWbPUATQGo6TqpNMhIACMDHcYAIAgBAYAIAiBAQAIQmAAAIIQGACAIAQGACAIgQEACEJgNB1LogMINDYT98zsBkmHJc24+91lt6cWWBIdQAaF3GGY2TEzu2Bmr/WV325mb5rZWTN7eFAd7v6Wux8YbUsbhiXRAWRQ1B3GCUnflfTDuMDMWpIek/RZSauSXjWzU5Jakh7p+/z97n6hmKY2CEuiA8igkMBw95fNrN1XfJuks+7+liSZ2TOS7nL3RyR9oYh2Nd7evdHXUEnlANCnzE7vayS93fN6daMskZnNmtnjkj5lZt8acNyimS2b2fLFixfza20dsSQ6gAzK7PS2hDJPO9jdL0k6uFOl7r4kaUmSFhYWUuuDNju2Dx+OvobauzcKCzq8ASQoMzBWJV3X8/paSe+U1JbmYkl0AIHK/ErqVUk3mtn1ZjYl6R5Jp0psDwBggKKG1T4t6RVJN5nZqpkdcPcPJT0k6UVJb0h61t1fz+l8+81saW1tLY/qAACSzL2+X/MvLCz48vJy2c0AgLFiZqfdfaG/nKVBAABBCAwAQJBaBgZ9GACQv1oGhrs/7+6LMzMzZTcFAGqjloEBAMgfgQEACEJgAACC1DIw6PQGgPzVMjDo9M4Z27gC0Bht0YqSsI0rgA21vMNAjtjGFcAGAgODsY0rgA21DAw6vXOUtl0r27gCjVPLwKDTO0ds4wpgQy0DAznqdKSlJWl+XjKLnpeW6PAGGohRUtgZ27gCEHcYAIBABAYAIEgtA4NRUgCQv1oGBqOkACB/tQwMAED+CAwAQBACA8NjNVugEZiHgeGwmi3QGNxhYDisZgs0BoGB4bCaLdAYtQwM5mEUiNVsgcaoZWAwD6NArGYLNEYtAwMFYjVboDEYJYXhsZot0AjcYQAAghAYAIAgBAYAIAiBAQAIQmAAAIIQGACAILUMDGZ6Vwgr2QK1UcvAYKZ3RcQr2a6sSO6bK9kSGsBYqmVgoCJYyRaoFQIDo8NKtkCtEBgYnUEr2dK3AYwdAgOjk7aS7ec/T98GMIYIDIxO2kq2L7xA3wYwhszdy27DyCwsLPjy8nLZzUC/iYnozqKfmbS+Xnx7AGxhZqfdfaG/nDsMFI9d+oCxRGCgeOzSB4wlAgPFY5c+YCyx4x7KwS59wNjhDgMAEITAAAAEqWVgsFotAOSvloHBarUAkL9aBgYAIH8EBgAgCIEBAAhCYAAAghAYAIAgBAaqg02VgEpjaRBUQ7cbbaIU75MRb6oksYQIUBHcYaAaDh9mUyWg4ggMVMP589nKARSOwEA1pG2eNDFBnwZQEQQGqiFpUyVJunIl2s417tMgNIDSEBiohv5NlVqt7cfQpwGUisBAdXQ60rlz0vp69EhCnwZQGgID1ZTWp5FWDmDkCAxUU1KfxvR0VA6gFAQGqqm/T2N+PnrNJD6gNMz0RnV1OgQEUCHcYQAAgoxNYJjZF83sB2b2YzP7XNntAYCmKSQwzOyYmV0ws9f6ym83szfN7KyZPTyoDnf/kbs/IOlrkr4ywuZiXLC6LVCoovowTkj6rqQfxgVm1pL0mKTPSlqV9KqZnZLUkvRI3+fvd/cLG3/+9sbn0GSsbgsUrpA7DHd/WdKv+opvk3TW3d9y98uSnpF0l7v/p7t/oe9xwSLfkfQTd/9FEe1GhbG6LVC4MvswrpH0ds/r1Y2yNN+U9BlJd5vZwbSDzGzRzJbNbPnixYv5tBTVw+q2QOHKHFZrCWWedrC7Pyrp0Z0qdfclSUuStLCwkFofxtzevdHXUEnlAEaizDuMVUnX9by+VtI7JbUF44aZ4EDhygyMVyXdaGbXm9mUpHskncqjYjPbb2ZLa2treVSHKolHRt17r/Txj0uzs8wEBwpS1LDapyW9IukmM1s1swPu/qGkhyS9KOkNSc+6++t5nM/dn3f3xZmZmTyqQ1XEI6NWVqI9Mi5dkn7zG+nkyWiVW8ICGClzr+/X/AsLC768vFx2M5CXdju532J+PgoMALkws9PuvtBfPjYzvQFGRgHlyhQYG30NlUcfRk2xRwZQqqx3GH8d/8HM/ijntuSGPoyaYmQUUKqsgfHixsS4jqQ7RtEgIFWWPTJ615mam4serDkFDCW409vMnpC0Jun3Jf27u//lCNuVCzq9G6bbjZYGWVmJAiXtv+3paYbgAgMM3ent7gckHZb0V5IumNn3c2wfMJzeIbdSelhIrDkF7NKOgWFmd5jZf5jZm5KelLTu7n/v7t8YffN2h07vBkpajHAQRlYBmYXcYRyV9OeS/kDRGk1/Z2ZfHWmrhkSndwNlDQBGVgGZhQTGf7v7v7r7/7j7TyX9iaKvpoDqyBIAjKwCdiUkMM6Z2d/0zMH4QNJ7I2wTkF3SkFvbWBB5dpY1p4AchASGS/qSpLfN7OeSzkp6ycxuHGnLgCyShtyePBl1fr/7bvRYX2fNKWAIWYbV/pak35V0S8/jBne/buAHS2Bm+yXt37dv3wNnzpwpuzkAMFbyGFb7f+6+7O5PuPufufsfVzEsJDq9EaB3Yh+T+YAgZe64B5QjnrMRD8NdWYleS3xdBQzAarVonqQ5G0zmA3ZEYKB5WCYd2BUCA83DMunArhAYaB6WSQd2pZaBwVpSGCjLMukAPsKe3gCALdjTG5CYfwEMgXkYaA7mXwBD4Q4DzcH8C2AoBAaag/kXwFAIDDQH8y+AodQyMBhWi0TMvwCGUsvAYLVaJGL+BTAURkmhWTodAgLYpVreYQAA8kdgAACCEBhAPPvbTNqzJ3pmFjiwDX0YaLb+2d9XrkTPzAIHtuEOA82WNPs7xixwYAsCA8220yxvZoEDHyEw0Gw7zfJmFjjwkVoGBjO9ESxp9ndsclL69a+3L4XOEuloKDZQArrdqK9iZUVqtaKO79lZ6b33pMuXN4+bnpbuu0968smt/R7T08wYR62kbaBEYABJ2u0oQPrFgdJvfl46d27UrQIKwY57QBZpnd1JYTHoeKBGCAwgSVpnd6uV7XigRggMIEnaUuiLiyyRjsYiMIAkaUuhHz3KEuloLAIDSNPpRB3Z6+ubHdrttnTvvdGfDx6Mnu+9l+G1aATWkgJC9K85tbIife97m++z9hQagDsMIMSgNadirD2FmiMwgJ10u8lzMpIwvBY1RmAAg8RfRYVieC1qjMAABgn5KirG8FrUHIEBDDLoK6ZDhxhei0apZWCwWi2G0rsa7UTKr8j8fDQn49w56eTJqIzhtai5WgaGuz/v7oszMzNlNwXjJu6zWFmR3JPXjur96qn/+Hh4LaGBGqplYAC7ltZn0Wolf/WUdDzDa1FTTNwDeqX1WayvR4/Q4xleixriDgPolTYsNq9yYIwRGECvtFVq04bLJh2ftrUrMOYIDKBX2iq1acNl+4+fnY2eL12iExy1wxatQJ7StnZlC1eMEbZoBYpAJzhqjMAA8kQnOGqMwADyEM8OX1mJ+jB6mUXlcQd470xyOsUxRpiHAQyrf3Ml9ygkep+lKDS+/vWo7PLlzTI2XsKY4A4DGFbSbG/3aHZ4/6CSDz7YDIsYM8MxJggMYFhpHdpJ61BlrQOoEAIDGFZah3arNXwdQIUQGMCw0maHLy4mzwKfmtp+LBsvYQwQGMCw0maHHz26vfz4cenYMTZewlhipjeQp2436sA+fz76munIEcIAYydtpjfDaoG89A+vZcgsaoavpIC8sJkSam5sAsPMPmlmj5vZc2Z2qOz2ANuwjhRqrpDAMLNjZnbBzF7rK7/dzN40s7Nm9vCgOtz9DXc/KOnLkrZ9twaUjnWkUHNF3WGckHR7b4GZtSQ9JukOSTdL+qqZ3Wxmv2dm/9z3+O2Nz9wp6eeSflZQu4FwWTdfAsZMIZ3e7v6ymbX7im+TdNbd35IkM3tG0l3u/oikL6TUc0rSKTP7F0n/OMImA9nFHduMkkJNldmHcY2kt3ter26UJTKzT5vZo2b2fUkvDDhu0cyWzWz54sWL+bUWCNHpRBslra9Hz4PCotuV5uai+RjxY26O1WtRWWUOq7WEstRJIe7+kqSXdqrU3ZckLUnRPIxdtg0YrW43Wrn2gw+2ll+6JN1/f/Rn7kxQMWXeYaxKuq7n9bWS3impLUCxDh/eHhaxy5cZiotKKjMwXpV0o5ldb2ZTku6RdKrE9gDF2WmoLUNxUUFFDat9WtIrkm4ys1UzO+DuH0p6SNKLkt6Q9Ky7v57T+fab2dLa2loe1QH522moLUNxUUGsJQWUIa0PQ4pWsz12jD4MlCZtLamxmekN1EqnE61cOzu7tXx2lrBAZbH4IFCWTodgwFip5R0GfRgAkL9aBoa7P+/uizMzM2U3BQBqo5aBAVRCtyu129LERDSDO57VvWdP9Nxuh83q7q0n6TPx+711x+eLP/Pgg4PrAAIwSgoYhf7NlNJMTw/eojWpnt7PhJ4n63nRaGmjpGoZGGa2X9L+ffv2PXDmzJmym4MmarejHfdCzM9H605lqSf+TJbzZDkvGq1RgRHjDgOlmZiQQn+3zKLFCrPUE38my3mynBeNxjwMoEhZZmoPOnanTZmGmRHObHJkRGAAo5C0mVKSnTZY2mlTptDzZD0vkIDAAEah04k6lefno69+Zmc3Z3W3WtHz/PzOHc/99fR/pvf93rrj88WfOXQovQ4gEH0YAIAtGtWHwUxvAMhfLQODmd4AkL9aBgYAIH8EBgAgCIEBAAhCYAAAgtQyMBglhUoatOpst7u5mm38uPrq6NFbNui9ubn0VWlDV7ztfT+0DI3BPAygCINWnZXS9/ce1vS0dN990pNPZlvxdnIyCqHLlweXseptLbH4IFCmQavOSrtfcTZEqyVduZJ87mFXvO2tB7WRFhjs6Q0U4fz5bOV5SgqL3nMP24Yi/g6ohFr2YQCVM2jV2VGvGhuvL5V07t7n3WLV28YgMIAiDFp19siRqH9gFKano/6JrCveTk5KU1M7l7HqbaMQGEARBq062+lIx49vrmYbu+qq6JEk6b3Z2eRVaY8eDV/xNn7/+HHp2LGdy+jwbhQ6vQEAW7BaLQBgKLUMDFarBYD81TIwAAD5IzAAAEEIDABAEAIDABCEwAAABCEwAABBCAwAQBACAwAQpJaBwUxvAMhfLQODmd4AkL9aBgYAIH8EBgAgCIEBAAhCYAAAghAYQNV0u9LcXLSrnVn052538HtJ5Q8+uLXsYx+L9veOX/c+ko5PKut9xHXt2RM9t9vR8e129HpiIvlzvZ/t/czERPQc/33ieuL64+err97+94g/l/bzDKm//73e4wfVFXrOYY4LlXd9/dy9to9bb73VgbHy1FPuk5Pu0tbH1JT7oUPJ77Va0aO/fFwfk5PR3zfr56ano59f/89zejq8/qT34nqT6go95zDHZflvJ6f6JC17wr+pbNEKVEm7La2sJL/XaklXrhTanLEzPy+dO7f5etDPM2u9UnJdoefc7XGhcqwvbYtWAgOokomJ6P8NsTtm0vr65uu8fp5m0XNSXaHn3O1xoXKsr1F7egNja+/e9PdareLaMa76f36Dfp5Z602rK/Scuz0uVN71JSAwgCo5ckSanNxePjUlLS4mv9dq1StMJiejv29W09PRz6/XkSNReWj9Se/F9SbVFXrOYY4LlXd9SZI6NuryoNMbY+mpp9xnZzc7LmdnNzsu095LKj90aGvZ1JT7xERyZ2/S8UllvY+4rrjDfX4+On5+Pnptlt5BHX+29zNm0XP894nrieuPn6+6avvfI/5c2s8zpP7+93qPH1RX6DmHOS5UTvWJTm8AQIhG9WGwWi0A5K+WgeGsVgsAuatlYAAA8kdgAACCEBgAgCAEBgAgCIEBAAhS63kYZnZRUv9qXDOS+sfbJpXNSXp3RE0bJKktRdQR8pmdjhn0ftp7IddjnK/FbuoJPX6312OYayGN9/XgdyPMvLt/Yltp0my+Oj8kLQWWJc50LKN9RdQR8pmdjhn0ftp7IddjnK/FbuoJPX6312OYazHu14PfjeEeTfxK6vnAsrLk0Zbd1BHymZ2OGfR+2ntVvh55tSNrPaHH7/Z6jOO1kPjdCG3LyNT6K6lhmNmyJ0yNR/G4FtXC9aiOoq9FE+8wQi2V3QB8hGtRLVyP6ij0WnCHAQAIwh0GACAIgQEACEJgAACCEBiBzOwqM3vSzH5gZp2y29NkZnaDmT1hZs+V3RZIZvbFjd+LH5vZ58puT5OZ2SfN7HEze87MDuVdf6MDw8yOmdkFM3utr/x2M3vTzM6a2cMbxV+S9Jy7PyDpzsIbW3NZroW7v+XuB8ppaTNkvB4/2vi9+Jqkr5TQ3FrLeC3ecPeDkr4sKffhto0ODEknJN3eW2BmLUmPSbpD0s2SvmpmN0u6VtLbG4ddKbCNTXFC4dcCo3dC2a/HtzfeR75OKMO1MLM7Jf1c0s/ybkijA8PdX5b0q77i2ySd3fi/2MuSnpF0l6RVRaEhNfznNgoZrwVGLMv1sMh3JP3E3X9RdFvrLuvvhrufcvc/lJT7V+f8w7fdNdq8k5CioLhG0j9J+lMz+56qtVxCnSVeCzObNbPHJX3KzL5VTtMaKe1345uSPiPpbjM7WEbDGijtd+PTZvaomX1f0gt5n3RP3hXWgCWUubv/r6SvF92Yhku7Fpck8Q9T8dKux6OSHi26MQ2Xdi1ekvTSqE7KHcZ2q5Ku63l9raR3SmpL03EtqoXrUR2lXAsCY7tXJd1oZteb2ZSkeySdKrlNTcW1qBauR3WUci0aHRhm9rSkVyTdZGarZnbA3T+U9JCkFyW9IelZd3+9zHY2AdeiWrge1VGla8HigwCAII2+wwAAhCMwAABBCAwAQBACAwAQhMAAAAQhMAAAQQgMAEAQAgMAEITAAApkZp8xs5NltwPYDQIDKNYtkn5ZdiOA3SAwgGLdIumXZvYxMzthZn9rZklLVQOVw34YQLFukXRB0aJx/+DuT5XcHiAYiw8CBTGzSUnvSlqR9A13f6XkJgGZ8JUUUJybFe1j8KGkKyW3BciMwACKc4ukf1O02c1xM/udktsDZEJgAMW5RdJr7v5fkv5C0rMbX1MBY4E+DABAEO4wAABBCAwAQBACAwAQhMAAAAQhMAAAQQgMAEAQAgMAEITAAAAE+X/22DgFA2zk2QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x345.6 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(6, 4.8))\n",
    "plt.plot(k, Pk, 'ro')\n",
    "plt.xlabel(\"$k$\")\n",
    "plt.ylabel(\"$p_k$\")\n",
    "plt.xscale(\"log\")\n",
    "plt.yscale(\"log\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "79af6940-9b8e-453a-8123-805c22b12d31",
   "metadata": {},
   "source": [
    "## 3、度保持的网络随机化"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "090e545a-62aa-4967-9ac0-93d4a5cd6122",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2018 2930\n"
     ]
    }
   ],
   "source": [
    "# 加载蛋白质网络数据集\n",
    "df = pd.read_csv(\"protein_interaction.csv\")\n",
    "G = nx.from_pandas_edgelist(df, 'source', 'target', create_using = nx.Graph())\n",
    "print(len(G.nodes()), len(G.edges()))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "71635af8-45c9-4a49-9e9d-b6efa8c44a69",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1647 2518\n"
     ]
    }
   ],
   "source": [
    "# 删除self-loops\n",
    "self_edges = []\n",
    "for e in G.edges():\n",
    "    if e[0]==e[1]:\n",
    "        self_edges.append(e)\n",
    "G.remove_edges_from(self_edges)\n",
    "\n",
    "Gcc = sorted(nx.connected_components(G), key=len, reverse=True) \n",
    "# 得到图G的最大连通组件\n",
    "LCC = G.subgraph(Gcc[0])\n",
    "\n",
    "# 获取原始网络的最大连通子图节点数和连边数\n",
    "N, M = len(LCC.nodes), len(LCC.edges())\n",
    "print(N, M)\n",
    "# nx.draw(LCC, pos=nx.kamada_kawai_layout(LCC), node_size=30, node_color=\"red\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "5426dd8f-d98b-4073-bb89-0dbc1a30aa88",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.05659957171711166\n",
      "0.015345050524785101\n"
     ]
    }
   ],
   "source": [
    "print(nx.average_clustering(LCC))\n",
    "\n",
    "# 生成度保持的随机化网络：调用double_edge_swap函数\n",
    "newG = LCC.copy()\n",
    "G_d = nx.double_edge_swap(newG, nswap=M, max_tries=5*M)\n",
    "print(nx.average_clustering(G_d))\n",
    "# nx.draw(G_d, pos=nx.kamada_kawai_layout(G_d), node_size=30, node_color=\"red\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "0eac5fc4-6cf0-4928-ad0a-af17d2c4a42c",
   "metadata": {},
   "outputs": [],
   "source": [
    "# 生成完全随机化的网络\n",
    "G_r = nx.gnm_random_graph(N, M)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "2f10b81e-b6c3-47a2-86e8-f8cd4eb21a05",
   "metadata": {},
   "source": [
    "### 测试小世界性质"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "8bf5ac10-d8b5-426f-b536-eeda6fd7fb69",
   "metadata": {},
   "outputs": [],
   "source": [
    "all_sp1, Pl = get_pdf_spl(newG)\n",
    "all_sp1_d, Pl_d = get_pdf_spl(G_d)\n",
    "all_sp1_r, Pl_r = get_pdf_spl(G_r)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "19eccaf3-5c98-4fac-97c3-c3ee7f21981d",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0, 0.5, '$p_d$')"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAF1CAYAAAD2uq3cAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABYq0lEQVR4nO3deXxU5fX48c8zWckCARKSkLBkkDUBspEQlgBat2q11VIXqnUraqt1q1a/LS7V+tNqrV3cqBVb69pal1atqCCLLJKwCGFnEiAbhD0hZJ3n98dkQgiTZBIyc+/MnPfrNa8kM/fOPYTMnLn3PM95lNYaIYQQoiMWowMQQghhbpIohBBCdEoShRBCiE5JohBCCNEpSRRCCCE6JYlCCCFEp4KNDsATYmNj9fDhw40OQwghfEZhYeEBrXWcq8f8MlEMHz6cgoICo8MQQgifoZTa3dFjculJCCFEpyRRCCGE6JQkCiGEEJ3yyxqF8C+NjY2UlpZSV1dndChC+Lzw8HCSk5MJCQlxex9JFML0SktLiY6OZvjw4SiljA5HCJ+ltebgwYOUlpaSkpLi9n5y6UmYXl1dHQMHDpQkIcQZUkoxcODAbp+dS6IQPkGShBC9oyevJUkUQvSib3/72xw5cqTTbR588EE+//zzHj3/l19+ycUXX9yjfdt79tlnqa2t7ZXnamvmzJkyj8nPSKIQ/uf112H4cLBYHF9ff93jh9RaY7fb+fjjj4mJiel021//+td861vf8nhMXfFEomhubu7V5xPmIIlCeJ4337hffx3mzoXdu0Frx9e5c8/4mM888wxpaWmkpaXx7LPPAlBSUsLYsWP5yU9+QmZmJnv37mX48OEcOHAAgEcffZQxY8Zw7rnnctVVV/H0008DcN111/Gvf/0LcHQReOihh8jMzGT8+PFs3boVgK+//popU6aQkZHBlClT2LZtW6fxvfrqq1x22WVccMEFjBw5kvvuu6/1sYULF5KXl0dmZiazZ8+mpqaGP/7xj5SXlzNr1ixmzZrFO++8w9133w3AH/7wB6xWKwC7du1i2rRpAHzxxRdkZGQwfvx4brjhBurr61v/Db/+9a+ZNm0a//znP1uPa7fb+dGPfsSvfvWrM/rdC+PJqCfhWc43bucnV+cbN8CcOd1/vjvvhPXrO3581SpoeQNrVVsLN94If/mL633S06Hlzd+VwsJCFixYwOrVq9Fak5uby4wZM+jfvz/btm1jwYIFPP/886fsU1BQwLvvvsu6detoamoiMzOTrKwsl88fGxvL2rVref7553n66ad5+eWXGTNmDEuXLiU4OJjPP/+c//u//+Pdd9/t+N8NrF+/nnXr1hEWFsbo0aO5/fbb6dOnD4899hiff/45kZGRPPnkkzzzzDM8+OCDPPPMMyxevJjY2FgqKyt56qmnAFi2bBkDBw6krKyM5cuXM336dOrq6rjuuuv44osvGDVqFNdeey0vvPACd955J+AYcrl8+XIAXnzxRZqampgzZw5paWn88pe/7DRuYX5yRiE865e/PJkknGprHfd7Qvsk0dX9bli+fDnf+973iIyMJCoqissuu4xly5YBMGzYMCZPnuxyn0svvZQ+ffoQHR3Nd77znQ6f/7LLLgMgKyuLkpISAI4ePcrs2bNJS0vjrrvuoqioqMs4zznnHPr160d4eDjjxo1j9+7drFq1is2bNzN16lTS09P529/+xu7dp7f0SUhIoKamhurqavbu3cvVV1/N0qVLWbZsGdOnT2fbtm2kpKQwatQoAH70ox+xdOnS1v2vuOKKU57v5ptvliThR+SMQnhUwu5V7CPhtPvjd1dS2ZMn7OSTP+C4tOXijZBhw+DLL3tyRLTWHT4WGRnZ7X3aCwsLAyAoKIimpiYA5s2bx6xZs3jvvfcoKSlh5syZbj9P2+fSWnPuuefy5ptvdrl/Xl4eCxYsYPTo0UyfPp1XXnmFlStX8rvf/Y7i4uJO923/e5gyZQqLFy/mnnvuITw8vMtjC3OTMwrhOZWVLpME0OH9Z+w3v4GIiFPvi4hw3N9D+fn5vP/++9TW1nL8+HHee+89pk+f3uk+06ZN4z//+Q91dXXU1NTw0UcfdeuYR48eJSkpCXDUH3pq8uTJfPXVV+zcuROA2tpatm/fDkB0dDTV1dWt2+bn5/P000+Tn59PRkYGixcvJiwsjH79+jFmzBhKSkpan+e1115jxowZHR73xhtv5Nvf/jazZ89uTX7Cd0miEJ6xZAlkZHj/uHPmwPz5jjMIpRxf58/vWT2kRWZmJtdddx05OTnk5uZy0003kdHFv23SpElccsklTJw4kcsuu4zs7Gz69evn9jHvu+8+HnjgAaZOnXpGI4ni4uJ49dVXueqqq5gwYQKTJ09uLZjPnTuXCy+8kFmzZgEwffp09u7dS35+PkFBQQwZMqS1kB0eHs6CBQuYPXs248ePx2KxcMstt3R67LvvvpvMzEyuueYa7HZ7j/8NwniqO6fIviI7O1vLOG6D2O3w1FPYH/glfxt0Lzfs+38dburun96WLVsYO3ZsLwXoPTU1NURFRVFbW0t+fj7z588nMzPT6LCEcPmaUkoVaq2zXW0vNQrRew4fhmuvZcl/j3FXzA7W7XO/l4w/mjt3Lps3b6auro4f/ehHkiSEz5JEIXpHQQG2797NfeV38i6XMSRa8+YLcNVVRgdmnDfeeMPoEIToFVKjEGdGa4498zL35y5mbNlnfBJ2KY8+Clu3Kq68EuLjXe/W0f1CCPPxSqJQSl2glNqmlNqplLrfxeOXKqW+UUqtV0oVKKWmubuvME7z0Rr+MvmvjLznOzxpv5errtTs2BXEr351cuBRZaWjFuGcI/d2zu/Q2nG/EMI3eDxRKKWCgOeAC4FxwFVKqXHtNvsCmKi1TgduAF7uxr7CAIsXlJAZX8rcr29i1LB61qy28+qb4Qwe7Hr70aPBQjNFu2RMvRC+xhtnFDnATq21TWvdALwFXNp2A611jT45/CoS0O7uK7xr50743qS9nH3DcI42RvDOvI0sLR5Kdk7nf0rh4XDWgEMUHUw4o1nSQgjv80aiSAL2tvm5tOW+UyilvqeU2gp8hOOswu19W/af23LZqqCqqqpXAhcnHTkCP7+riXGjm/i8IIbHh73E1h3BzP71eNxtb59qrWMTqbBjh0dj9bSHH364tcGfcO2mm25i8+bNRofhUtvGjUbypTbv3kgUrt5GThtBr7V+T2s9Bvgu8Gh39m3Zf77WOltrnR0XF9fTWEU7TU3w4oswckQzzzxr4Rr739h+67M8sOMGwq0dXGfqQFpGCDs5i/oNWz0ULSQkOObZtb8leGgi+JnwVktuTx2ns+d9+eWXGTeud68SO1u5+wtfavPujURRCgxp83MyUN7RxlrrpcAIpVRsd/cVPdfRG2x4ONx6K4w7torCqJn89f1YEp+fB91YmN0pdfoAmglm23LPnfHt29e9+931m9/8htGjR/Otb33rlJbfu3bt4oILLiArK4vp06e3znretWsXkydPZtKkSTz44INERUUBjoWHZs2axdVXX8348eNpbm7m3nvvZdKkSUyYMIGXXnqp9bmfeuqp1vsfeughl3E9/PDDXHPNNZx99tmMHDmSv7R0yHX3OBUVFeTn55Oenk5aWlprs0NXrcnh1Jbiv/3tb8nJyWmNpaSkhAkTJgCnfqqNioril7/8JRMnTmTy5Mnsa/nP6Oh31JarVu633nor2dnZpKamnvJ76ahl+8GDBznvvPPIyMjg5ptvPqUPV0ft48eMGcNNN91EWloac+bM4fPPP2fq1KmMHDmSr7/++rQ4/b7Nu9baozccczVsQAoQCmwAUtttcxYnZ4lnAmU4zia63NfVLSsrS4vucYxNcn17l+9pe3qG1rt2ndExNm50PN8bOb/v1n6bN29u/f6OO7SeMaPjW2f/jo72ueOOzo9fUFCg09LS9PHjx/XRo0f1iBEj9FNPPaW11vrss8/W27dv11prvWrVKj1r1iyttdYXXXSRfuONN7TWWr/wwgs6MjJSa6314sWLdUREhLbZbFprrV966SX96KOPaq21rqur01lZWdpms+lPP/1U//jHP9Z2u103Nzfriy66SC9ZsuS02B566CE9YcIEXVtbq6uqqnRycrIuKytz+zhPP/20fuyxx7TWWjc1Neljx47pqqoqPX36dF1TU6O11vqJJ57QjzzyiNZa62HDhuknn3yy9fgTJ07Uu1r+Lp544onWY8yYMUOvWbNGa601oD/88EOttdb33ntv6zYd/Y7aKi4u1kopvXLlytb7Dh482BrvjBkz9IYNG1pj++Mf/6i11vq5557TN954o9Za69tvv701/v/+978a0FVVVa3/rzU1Nbq6ulqPGzdOr127VhcXF+ugoCD9zTff6ObmZp2Zmamvv/56bbfb9fvvv68vvfTS0+JcsGCBTklJ0UeOHNEnTpzQQ4cO1Xv27Onyd1lVVaW11rqiokJnZ2drrbW+/PLLdXZ2ti4tLdWvvvqqvv/++/WJEyd0cnKy3rZtm9Za62uuuUb//ve/d/l/MmPGDL1y5Up95ZVXtv7fttf2NeUEFOgO3lM9fkahtW4CbgM+BbYA72iti5RStyilnM1iLgc2KaXW4xjldEVL7C739XTM4lSXzY1DrVwBLZ9yemrUKAhWTWyy9emlyLxj2bJlfO973yMiIoK+fftyySWXAI4WHStWrGD27Nmkp6dz8803U1FRAcDKlSuZPXs2AFdfffUpz5eTk0NKimPW+sKFC/n73/9Oeno6ubm5HDx4kB07drBw4UIWLlxIRkYGmZmZbN26lR0d1Hac7cxjY2OZNWtW6yded44zadIkFixYwMMPP8zGjRuJjo7usjV525biP/jBD3jnnXcAePvtt09rNw4QGhraunxr21bqnf2O2mrfyv2dd94hMzOTjIwMioqKTqmFuGrZvnTpUn74wx8CcNFFF9G/f3+g8/bxKSkprT2tUlNTOeecc1BKMX78+Nbnbc+f27x7ZWa21vpj4ON2973Y5vsngSfd3Vd4WZvLIWciNBRGDTzoGPnU2Nijy1dddRnvrLDewy7jLc97+hPb7XZiYmJY39lCSi60bcmtteZPf/oT559//inbfPrppzzwwAPcfPPNp9z/3HPPtV5e+vjjj13G5vzZneOA4430o48+4pprruHee++lf//+nbYmb/u8V1xxBbNnz+ayyy5DKcXIkSNP2z4kJKQ1prat1N3V9njFxcU8/fTTrFmzhv79+3PddddRV1fX+rirlu3g+v9Pd9JsrG3LdovF0vqzxWLpMH5/bvMuM7OFV6Va6yjS4xzjbH1Efn4+7733HidOnKC6upr//Oc/APTt25eUlJTW68JaazZs2AA42ns7V6R76623Onzu888/nxdeeIHGxkYAtm/fzvHjxzn//PN55ZVXWmsDZWVl7N+/n5/+9KesX7+e9evXM7hl0soHH3xAXV0dBw8e5Msvv2TSpEluH2f37t0MGjSIH//4x9x4442sXbu209bk7Y0YMYKgoCAeffRRl2cTnXH3d9TWsWPHiIyMpF+/fuzbt49PPvmky33y8/N5vWUp3E8++YTDhw+33t/d9vHd5S9t3iVRCK9KzQhhFyM4sXaLR57fEy1DMjMzueKKK0hPT+fyyy8/5c3k9ddf569//SsTJ04kNTWVDz74AHCMaHnmmWfIycmhoqKiwxbjN910E+PGjSMzM5O0tDRuvvlmmpqaOO+887j66qvJy8tj/PjxfP/73z/lTaWtnJwcLrroIiZPnsy8efNaE4g7x/nyyy9JT08nIyODd999lzvuuKPT1uSuXHHFFfzjH//gBz/4QXd+rW7/jtqaOHEiGRkZpKamcsMNNzB16tQu93nooYdYunQpmZmZLFy4kKFDhwI9ax/fXX7T5r2j4oUv36SY3X3xlv0uC8Dxlv29epx//qNOg9Zrb3nJ7X1cFd7M7vjx49put2uttX7zzTf1JZdc4pHjPPTQQ62FdV/jrd+ROF13i9nSPVYAUPn3hbx9/f+4svE1NpJGGkWOhk3z5wM9X/SnvbQsx3XcTesaMWBZI68pLCzktttuQ2tNTEwMr7zyitEhmY78jnyHJArhMGcOtgUJ8AUMZ7djZbjf/OaMVoZz5ayzINTSSJHNv3s+TZ8+vbVe4UkPP/ywx4/hKd76HYkzJ4lCtLLZhzOIfURtXAlpaR45RnAwjB54gKIDCY5p38HyJyiE2UkxW7Sy7bZgxXbG8yW64hj5NBZ27fLocYQQvUMShWhVvD8Ka3jFycUkPCQtM5RirNQUbut6YyGE4SRRCMAx/21PzQCsscc8fqzU6QMA2LLM+A6eQoiuSaIQAOzdC80EYR165pNzupKa7WjhUbSuwePHMitpdS18iVQSBQC2rQ1AKNZRnv+TsFoh3FJPkQd6PmW8lMH6yvWn3Z+ekM66m9ed8fM7x5VbLP7xGevZZ5/lhz/8IRG9eLnRW+3Thff4x1+7OGO2goMAWNP7evxYQUEwJvYAmw4kQC+/qeQl5xEaFHrKfaFBoUxJntLj55RW1yZpdS0MI2cUAgDbxuOEUs/gTO+s8JM24gRL9o+F4mLH5Ao33fm/O12eMTjVN9XTZD/18lmTvYl1leuY+epMl/ukJ6Tz7AXPdnrcbdu2sWDBAp5//nnAsT7FgAEDaG5u5pxzzuGbb75pXYshNjaWtWvX8vzzz/P000/z8ssv88gjjzBt2jQefPBBPvroI+bPnw84Jp0tWLCA1atXo7UmNzeXGTNm0L9/f3bu3Mk///lP5s+fz6RJk3jjjTdYvnw5H374IY8//jjvv//+aXGuX7+edevWERYWxujRo7n99tvp06cPjz32GJ9//jmRkZE8+eSTPPPMMzz44IM888wzLF68mNjYWCorK3nqqacAR8fcgQMHUlZWxvLly5k+fTp1dXVcd911fPHFF4waNYprr72WF154gTvvvBNwtJlYvnw5AC+++CJNTU3MmTOnV7uYCmPIGYUAwLZTM5wSgkaN8MrxUjNC2ctQjq3p3ZFPYcFhxEfGo1oWR1QoEiITTjvL6C5pdW18q2thHDmjEADYysNICdoGg0Z55Xip+bHwPGxeeoDJV7m/X1ef/AEqqiuw/tFKXVMd4cHhFN5cSELUmZ0pSatr41tdC+PIGYUAwHa4P9Z+hzpf0KEXpU5yFE83rWvs9edOjE7k+vTrsSgL16dff8ZJoj1pdW1Mq2thHEkUgiNH4HBjNNbEE1475vDhEBFUR1GxZ1a7m5c/j2lDpzFvxrxef25pdW1Qq2thGNXZqa+vys7O1jKG233rCu1kZlt495JXueyD67x23EkJe+hftZ2FjWdDJ8NNt2zZwtixYz0ez+aqzdQ2nj6nICIkgnFx4zx+fCG8xdVrSilVqLXOdrW9nFEIbIWOyyDWcd5dyzrVWkeRfSx0UJj1tsiQyNYiuJNCERkS2cEeQgQGSRQC2/qjAKRMivXqcVMzwygnicNf7/DqcTuSGJ142n1KKQZHn75inBCBRBKFwLa1gYEcoN+EYV49btqMgQAULTG+lcXxhuPYDtvQnLwUq1AM7DOQkKAQAyMTwniSKAS23cGO9uItBVZvSc2NAqBofdc9nzxVS2tsbqTkSAlbDmyhrqmO5L7Jp1x+krMJ4W968lqSeRQC2/4osiK2QOiZTUrrriFDIDroeJc9n8LDwzl48CADBw50ORehJ+zazv7j+ymvLkdrTXxkPInRiQRbgqlvqqeqtoqw4DA5mxB+RWvNwYMHuz2vRRJFgGtuhpLjscweWt31xr1MKRgXV8Wm/Ylgt3c48ik5OZnS0lKqqqp65bgnGk9w6MQhmuxN9AnpQ//w/tQcrWEHjlpJs72Zg9UHCbGEsOXgll45phBmER4eTnJycrf2kUQR4EpLoUkHk+KF9uKupI2o4z+VYxx9zoe5rpGEhISQkpJyxsfaemArd396N5/s/ITRA0fz+/N/z4yRrieM/f2zv/OH1X+g+oHqM27/IYSvkxpFgLNtcswbsI425s0wNSOU/cRzYMV2jx3jSN0R7vn0Hsa/MJ6v9n7F7877Hd/c+g0Xjryww32yErNoaG6gaH+Rx+ISwldIoghwxWscI46s6dGGHD91pmNIbtHSg73+3M32Zv5S+BdG/WkUv1/1e66beB07bt/B3Xl3d3mWkDU4C4DCisJej0sIXyOJIsDZNtUSRBNDJp0+h8AbUic71r/YtL53L30t272MSX+ZxNz/zmV07GgK5hbwl0v+wqDIQW7tP6L/CPqF9aOwXBKFEFKjCHC2XXaGsZvg0d5pL97e4MEQE1zda6vd7T26l/s+v4+3Nr1Fct9k3rz8Ta5IvaLbo6WUUmQmZsoZhRBIogh4tvJwrMF7oJ8xiUIpSI3bT9H+eNDare61HS13mhCZwNH6o2g0D+Y/yH1T7yMytOftN7ISs/jT13+isblRhsmKgCaXngKc7fAArP0PGRpD6oh6iprHoPeWurW9q+VOASqPV3LxqIvZ8tMtPDLrkTNKEuCoU9Q311NUJQVtEdgkUQSw6mqoaozBOriu6409KDUjlIPEsu+rnW5tPy9/HhZ16p+uQvHuD97lndnvMDxmeK/ElZXYUtCWOoUIcJIoAljxDkcB2Wr1zmJFHUmbFQe4P/LJuTBRsMVx5TRIBXFz1s1cNvayXo1rxIAR9A3rK3UKEfAkUQQw29ctQ2PTIgyNI3VKPwCK1ru/2t28/HmtPWtCg0J5aOZDvR6XRVmkoC0EkigCmm1dS3vx7IGGxjFoEAwMPtqtkU+J0YmtQ109sdypU1ZiFhsqN9DY3PtLtgrhKyRRBDDbtkb6cYT+6d5tL96eUpA6aD+bDiY4Rj65wa7tHKs/RmJUokeWO3XKSnQUtDdXbfbYMYQwO0kUAax4bxBWVYxKTjI6FNJG1DlGPpVXuLX9zkM7Od54nMfOfsxjZxNwcoZ2QbksrSsClySKAGarisYaub/T9aq9JTUzlKPEUL7UvZFPzjdu58gkTzlrwFlEh0ZLnUIENK+8QyilLlBKbVNK7VRK3e/i8TlKqW9abiuUUhPbPFailNqolFqvlJKPdb3Ebofimjisg7zfXtyV1JmOeoO7I58KywsJCwpjXNw4T4YlBW0h8EKiUEoFAc8BFwLjgKuUUu1f3cXADK31BOBRYH67x2dprdO11tmejjdQVJRr6nUY1qHNRocCQOrUGMD9nk+FFYVMTJjolRnTUtAWgc4bZxQ5wE6ttU1r3QC8BVzadgOt9Qqt9eGWH1cB3VtVQ3Sbc8STdbQ5WlPExiniQw5SVNz1yCe7trO2Yi3Zid753OCcoS0FbRGovJEokoC9bX4ubbmvIzcCn7T5WQMLlVKFSqm5HogvINnWOC7xWDP6GRzJSalxVRQd6Hrk085DO6luqG4tNHta6wxtufwkApQ3EoWrab8u3wmUUrNwJIpftLl7qtY6E8elq58qpfI72HeuUqpAKVXQW0tm+jPbplosNDM015j24q6knlVPUfNodOW+TrfzViHbaeTAkY6CtrTyEAHKG4miFBjS5udkoLz9RkqpCcDLwKVa69aKpta6vOXrfuA9HJeyTqO1nq+1ztZaZ8fFxfVi+P7JZoMh7CV09JkvMdpbUjPDqCGaPV/aOt3OW4VsJ4uykJGYIWcUImB5I1GsAUYqpVKUUqHAlcCHbTdQSg0F/g1co7Xe3ub+SKVUtPN74Dxgkxdi9nu2inBSQsuhT++sA9Eb0txc7c6bhWynrMQsNuzbQJPdmLXFhTCSxxOF1roJuA34FNgCvKO1LlJK3aKUuqVlsweBgcDz7YbBxgPLlVIbgK+Bj7TW//N0zIHADO3F20vNd7QS6aznk7cL2U7Zg7Opa6qTgrYISF5ZuEhr/THwcbv7Xmzz/U3ATS72swET298vzkxtLVQ2DsQ6uN7oUE4R018xOLSq05FP3i5kO7VtOT4hfoJXjy2E0Yyfkiu8rmSrY/0J6whj24u7khpXxaaDHRfYvV3IdmotaEudQgQgSRQByLZqP2B8e3FX0s6qY0vTSOz7XI9c83Yh20kK2iKQSaIIQLb1xwCw5sQaHMnpUjPDqCWSki92uXzciEK2k3OGthS0RaCRRBGAbNsaiaKa2Cxj24u7kjqzZbW7ZaePfDKqkO2UlZjFiaYTbKnaYsjxhTCKJIoAZNsbjNVSgooz3xnFuBmORLFp/ek9qIwqZDs5jyuXn0SgkUQRgGxVfbFGVTlWDDKZvv0UQ0MrXI58MqqQ7TRq4CiiQqNkhrYIOJIoAozWYDs+COugGqND6VDqoAMUHTx9MSKjCtlOFmUhI0EK2iLwSKIIMPsq7JzQfUgxSXtxV1LPqmdL00ia9x045X4jC9lOWYlZrK9cLwVtEVAkUQQY29eON1/r2DCDI+lYamYY9YSz64uS1vuMLmQ7ZQ12FLS3HthqaBxCeJMkigBTXNDSXjy9r8GRdCxtVsvIpzY9n4wuZDu1naEtRKCQRBFgbEUnABg+ZbDBkXRs7Mx4AIo2nOz5ZHQh22nUwFFEhkS2xiNEIJBEEWBsNkiilPCRQ7re2CCRUYqUsDKKik/OHDe6kO0UZAmSGdoi4EiiCDC2ij5Yw8ohxBxLoHakfc8nMxSynaSgLQKNJIoAYzs6AOuAI0aH0aW0kXVsaxpB4/7DpilkOzlnaEtBWwQKSRQBpK4OyhrisCbVGR1Kl1Izw2kklJ2fFZumkO3UOkNbCtoiQEiiCCC7N1WjsWAdYf7/9tRZgwAoWnbINIVsp9EDRxMZEil1ChEwvLJwkTAHR3vxaKzjI40OpUtjZiZgoZlN65uoNkkh2ynIEkR6QrokChEwzP/RUvQaZ3vxlJw4gyPpWp9IC9awMopKIkxVyHZyFrSb7ead4S5Eb5FEEUCKtzcSzgkScoYaHYpb0gZVsengIFMVsp2yBmdR21grBW0RECRRBBDb3hCsQbtR/cw7K7ut1LPq2dGv2VSFbKfWGdpy+UkEAEkUAcR2IBprtOslRs0oNTOM5sEbAPMUsp3GxI4hIiRCRj6JgCCJIkA42ovHYx103OhQ3JY6axAkFhKiQ0xTyHaSgrYIJJIoAsTBykaqdTTW4b5TfB09azAMLiDuqNVUhWyn7MRs1lWuk4K28HuSKAKEbeU+AKxjzNtevL2QPgpL4lrCSs11NuHkLGhvO7jN6FCE8ChJFAHCVnAIAGtmjLGBdMPOQzuxh9VwvDjX6FBckpbjIlBIoggQzvbiKVPN2168PeeM7KrSc6mrqjY4mtO1FrSlTiH8nCSKAGErhgQqiLCevha1WRWWFxKiQ9FV49n66W6jwzmNFLRFoJBEESBsFRGk9KkEi+/8lxdWFDImaizYQyhadsjocFzKSsxiXYUUtIV/8513DXFGin2kvbiTs7X41FFTCKGBTRvMufZDVmIWxxuPs/3gdqNDEcJjJFEEgMYGzZ7GBKxJ9UaH4jZna/FJQ7IZFb6HomJzNjJsbTkul5+EH5NEEQD2rD+EnSCsZwUZHYrb2rYWTxtURdEhc9ZWxsSOoU9wHxn5JPyaJIoA0DqHYkKUwZG4r+0a2aln1VPcNITj+803qzzYEiwFbeH3JFEEANsGx9BSa06swZG4r21r8dSscDQWtny6x+iwXMpKzGJtxVopaAu/JYkiANi2NxFKPYNzhxgdilvar5HddrU7M8oaLAVt4d8kUQQAW2koKcF7sUSEGx2KW9qvkT1i1lDCqKPIxCOfQArawn9JoggAtgN9sUYfMDoMt7VfIzs4PJgx4SUUlUQYGVaHxsaNlYK28GuSKAKArTYea7z5CsEdKXSxRnbqoAMUHUo0MKqOBVuCmZgwUc4ohN+SROHnDpfVckTHkDLcbnQobnO1RnbqyAZ2NyVTva/WwMg6lpWYxbrKddi17/yehXCXJAo/V7yiAgDrWN+oT7QvZDulZTnao2/+dK8RYXUpKzGLmoYaKWgLv+SVRKGUukAptU0ptVMpdb+Lx+copb5pua1QSk10d1/ROV9rL96+kO3UOvJpuXlHPoG0HBf+yeOJQikVBDwHXAiMA65SSrVfiaYYmKG1ngA8Cszvxr6iE7bNdYDvtBd3vtG2XyM7ZeYw+lBL0QZzzlUYFzeO8OBwqVMIv+SNM4ocYKfW2qa1bgDeAi5tu4HWeoXW+nDLj6uAZHf3FZ2zFSti1QH6Dh9gdChuKSgvOK2QDWAJD2VseDGbTNrzKdgSzMR4KWgL/+SNRJEEtL2wXNpyX0duBD7p4b6iHVtlBNY+FaCU0aG4pbCikPSEdJdrZKfFV1F02Jwjn+Bky3EpaAt/441E4eodSrvcUKlZOBLFL3qw71ylVIFSqqCqqqpHgfoj27GBWAccNToMtzgL2e0vOzmlntVAWVMCRyrrvByZe7IGZ1HdUM2OgzuMDkWIXuWNRFEKtO0dkQyUt99IKTUBeBm4VGt9sDv7Amit52uts7XW2XFxcb0SuK9rqm9md+NgrMkNRofilo4K2U6pWY6RW5sXlnozLLfJDG3hr7yRKNYAI5VSKUqpUOBK4MO2GyilhgL/Bq7RWm/vzr6iY6WF+2giBOtI3xgF3VEh2yn17HgANi0158in1oK2jHwSfibY0wfQWjcppW4DPgWCgFe01kVKqVtaHn8ReBAYCDyvHNfSm1rODlzu6+mY/YVt1X5gMNYJ0UaH4paOCtlOQ/OHE0U1Rd+Yc+RTSFCIFLSFX/J4ogDQWn8MfNzuvhfbfH8TcJO7+wr3FH/jaC+ekjvI4Ejc01khG8DSJ4xx4VspKjHnyCdwnA299s1r2LUdi/KNMzkhuiJ/yX7MtqOZYBpJnmTekUJOXRWynczc8wlOFrR3HtppdChC9BpJFH7MVhrCsOBygsO9cuJ4RroqZDuljmqgsjmOg+XmXP+7taAtdQrhRyRR+DHbwX5Y+/lGe/GuCtlOaS0jn4oWlnk8pp4YFzeOsKAwqVMIvyKJwo/ZahN8pr14V4VsJ+fIp6LlhzvdzighQSHSclz4HUkUfurYniMc0LFYh7ucn2g6XRWynZKmW+nLUdOOfIKTa2jLDG3hLyRR+KnW9uLjzN9e3N1CNoDqE05q+C7T9nwCR6I4Vn+MXYd2GR2KEL1CEoWfshW2tBfP6m9wJF1zt5DtlBZ/wNw9n1r+Hc4lXYXwdZIo/JRts2NUkHWqed9QndwtZDuljmzgQPMA9pc1ejKsHkuNS5WCtvArkij8lK1E0V8dJmaI+Wdlu1vIdkrN7gNA0Wcu234ZLiQohAnxEyRRCL8hicJPFe+LICVin9FhuMXdQraTc7W7TcvM2fMJpKAt/IskCj9lOxqLdaD524t3p5DtlDB1BAM4aO6RT4OloC38hyQKP2Sva6C4KRlrkvnbi3e3kA2gIiNIDd9FUUmUByM7M9JyXPgTSRR+qHxNGQ2EYR1l/tYd3S1kO6XGH6TocCLapNNEUgelEhoUKq08hF+QROGHHO3FwTrR/wrZTqkjGzjc3I+KvU0eiuzMhAaFSkFb+A1JFH7I9k0NANZc86/0191CtlNadkvPp88rPBFWr3AWtLVZT3uEcJMkCj9k29GMhWaGZpt7HYqeFLKdzN7zCRyJ4mj9UXYdloK28G2SKPyQrSyMoSEVhIQqo0PpVE8K2U5xU0YSx37Tj3wCaTkufJ8kCj9kO9QPa7+DRofRpZ4WshMSQEVFUsUgXi7MQClQynG/maQNSnMUtKVOIXycJAp/ozW2E4lYE2qNjqRLPS1k7+tgHmFH9xslNCiU8YPGS6IQPk8ShZ+pLdnPPh1PynCjI+laTwvZvkQK2sIfSKLwM8VfOfofmb29+JkUsn1J9uBsjtQdwXbYZnQoQvSYJAo/Y1t7BABr9gBjA+nCmRSyfUlrQVsuPwkfJonCz9i2tLQXn2Kyym47PS1k+5rWgraMfBI+TBKFn7GVKKJVNQMHhxkdSqd6WsgGiO93olv3G0kK2sIfSKLwM7Z9UVgj9qHMPYXijArZlUf6oP/xOnrYcH7PnQBU/PldKo/06eUoe4cUtIWvk0ThZ2zVsVhjzd1evFcK2XPmQEkJWffMAqCwz7Reiq73ZQ3O4nDdYYqPFBsdihA9IonCj+ia49iahmIdYs4lQp16s5CdcckQFHYKPz3QC5F5RmvLcalTCB8licKPVK7ZSx19sI4y97yE3ixkR+WmMoZtFK4177W2tEFphFhCpE4hfJYkCj9iW10FgHWieRf0gTMrZJ8mLIysATYK9pq3AWJYcBjj48dTUF5gdChC9IgkCj9SvNHRXjwlN97gSDrX2zOys0bVUF4fS2WFeYvFUtAWvkwShR+x7WxGYWfYhH5Gh9IhT8zIzprmGO1U+LHJmj21kZUoBW3huyRR+BFbaRhJIfsJ72Pe6/WemJHdWtBeaM6OuRkvZXDLR7cAMOKPI1CPKNQjioyXMgyOTAj3SKLwI7ZDMVj7HTI6jE55YkZ2VM44Uxe085LzCA0KPeW+0KBQpiRPMSgiIbpHEoW/aG7GVmf+9uK9Wsh2Cgsja0CxaQva8/LnYVGnvtSCVBDzZswzKCIhusftRKGUmqOUGuXJYETP1e0spYxkrFajI+mcp1qLZ402b0E7MTqR69Ovb00WoUGhXJ9+PQlR5u7HJYRTd84oqoDnlVJLlVL/Vko94amgRPeVrKoEwJpqzjYW4NnW4llTHW3VCz+u7PXn7g3z8ucRpIIAUCg5mxA+xe1EobVeCKzWWucDPwLMPVg/wNjWHgbM3V7ck63FMy4dauqCdmJ0IpePvRyAKUOmyNmE8CndrVH0VUplAfVApAfiET1k29IAgHWyOa/Tg2dbi58saJu37Pb0eU8TpIJI7ptsdChCdEt3X1V3A1OBl4D/9X44oqeKd1uIULUMSgwyOpQOeaSQ7RQaauqCNkBS3yTOTjmbb/Z9Y3QoQnRLl4lCKXWhUmq1Umob8DqwRmt9vdb6bXcPopS6QCm1TSm1Uyl1v4vHxyilViql6pVSP2/3WIlSaqNSar1SSnogdMC2P5KUyP2mbi/u6TWyzVzQdspNymXj/o0cbzhudChCuM2dM4rncZxJTAbmA08ppa5y9wBKqSDgOeBCYBxwlVKq/UfKQ8DPgKc7eJpZWut0rXW2u8cNNLZjcVhjjxkdRoe8sUb2yRna5ixoA0xOnoxd26VBoPAp7iSKfVrrr7TWh7XWnwPnA7/sxjFygJ1aa5vWugF4C7i07QZa6/1a6zWAuftjm5Q+dBibfZip24t7Y43sky3HzVnQBshJygFgVekqgyMRwn3uJIoSpdRjSinn1NJGoLobx0gC9rb5ubTlPndpYKFSqlApNbcb+wWMA4W7qSHa1O3FvbFGdmtBe515r7/FRcZh7W9lddlqo0MRwm3BbmyjgcuAHyuldgBDgdeVUiO11jvc2N/Vq7Y7F5Gnaq3LlVKDgM+UUlu11ktPO4gjicwFGDp0aDee3ve1thdP72twJKfLeCmD9ZXrW39Ofynd8TUhnXU3r+vdg7UUtBftzend5+1lk5Mn82XJl0aHIYTbujyj0FpfpbUeBwwD7gQewTE09mWl1N7O9m1RCgxp83MyUO5ugFrr8pav+4H3cFzKcrXdfK11ttY6Oy4uzt2n9wu2TY62HdbJ5vt3e7vPka8UtMuryyk9Vmp0KEK4pTsT7uq01gVa679qrX+mtZ6htR7S9Z6sAUYqpVJaLl9dCXzozjGVUpFKqWjn98B5wCZ3Yw4Utp12AIaPM9/UFm/3OWotaH9k3oJ2blIuAKtL5fKT8A0en52ktW4CbgM+BbYA72iti5RStyilbgFQSiUopUpxjK76lVKqVCnVF4gHliulNgBfAx9prWX+Rju28jASQw4QEWF0JKdLjE7k4pEXt/7s6T5HZp+hDY7LbqFBoVLQFj7DnRrFGdNafwx83O6+F9t8X4njklR7x4CJno3O99kO9ccacxCINToUl8bGjXV8RMDzXVOjJo01fUE7LDiMjIQMKWgLn2HefgfCPfX1FNcnYk2sMzqSDhVWFBITFoNFWTzfNdU5Q3uPuZeDzU3KpaC8gCZ7k9GhCNElSRQ+rmHHbvYyhBSTthdvbG5k6e6lXDL6EqYNneaVrqlZY45T3mDugvbk5MmcaDrBxn0bjQ5FiC5JovBxe1ZXYCfItO3FCysKqWmo4Tujv8OS65Z4pWtqa8vx/1Z4/Fg9lZvcUtCWy0/CB0ii8HG2tUcAsGYPNDaQDiwqXgTAzOEzvXbM1oL2Z+YtaKfEpBAbESuJQvgESRQ+zra1pb24SdehWFyymAnxE4iN8F6hPWrSWMaobRSuM++ft1KKycmTZeST8AnmfSUJt9h2BxGm6kkcbL5RPvVN9Szfs5xZw2d598A+VNDeemArR+qOGB2KEJ2SROHjbFVRpETux2LC/8nVZaupa6rj7JSzvX7srNHmL2g7J96tKVtjcCRCdM6Eby/CbVpjqx6ENa47PRq9Z1HxIizKQv6wfK8fu3WGtokL2jlJOSiUXH4SpieJwofpikp26RSsQ8w5Fn9xyWIyEzOJCY/x+rF9YYZ2v/B+jIkdIwVtYXqSKHzY4Q17OEY/rGPM1168trGWlXtXer8+0cIXCtrgmE+xumw1Wpv3EpkQ5n4ViU4Vf93SXnyi+dqLr9i7gkZ7oyH1CQBCQsgaUELBXvMXtA/UHsB22GZ0KEJ0SBKFD7Ntcqy7nJI7yOBITreoeBHBlmCmDZ1mWAxZo2scBe1yu2ExdEUm3glfIInCh9l2OS5XpJhwZbvFJYvJScohKjTKsBh8oaCdNiiNiJAIaTkuTE0ShQ+zlYUTF3KY6GijIzlVdX01a8rWGFafcMr47rCWGdqHDI2jM8GWYLIHZ7OqTEY+CfOSROHDbIf7Y405bHQYp1m2ZxnNutm4+kSLqOwxLQVt801GbCs3KZf1leupb6o3OhQhXJJE4atqarA1JmNNPGF0JKdZVLyI0KBQ8pLzjA3ERwrak5Mn09DcwLrKXl5DXIheIonCF73+Ok0jx7KbYVh3fAqvv250RKdYXLKYKUOm0CfE+I62jpbjceYuaMvSqMLkJFH4mISYE6gfziGkci/NBPObE3ejfjiHhBhznFkcOnGIdRXrDK9POPlCQTupbxJJ0Uky8kmYliQKH7PvqOtP6R3d721Ldy9Fow2vTzj5QstxQDrJClOTRCF61aLiRUSERJCTlGN0KEDbgra5/9Rzk3IpPlJM1fEqo0MR4jTmfvUIn7O4ZDHThk4jNCjU6FAcfKSgLRPvhJlJohC9Zv/x/Wzav8k09QknXyhoZyVmEaSC5PKTMCVJFKLXfFnyJYBp6hNOWdMjAHMXtCNDIxkfP17OKIQpSaLwMfEdXEHp6H5vWlS8iOjQaDITM40O5RS+0HIcHHWKr8u+xq7Ne+YjApMkCh9TWQlPDHsBgKoq0Npxq6w0ODAc9YkZw2cQbAk2OpRTtBa015v7z31y8mSO1R9j64GtRocixCnM/coRp2toYMXeZEb1309srNHBnFR2rIztB7ebrj4BQHAwWQN9oKAtE++ESUmi8DF6/QZW2nOZMvG40aGcYnHJYsB89QknXyhoj44dTb+wflKnEKYjicLH7PrvFqoYRN4FMUaHcopFxYsY0GcAE+InGB2KS1nTWgra/yk3OJKOWZSFnKQcGfkkTEcShY9Z+VkNAFMu6m9wJKdaXLKYmcNnYlHm/JPyhZbj4Lj8tHH/Ro43mOuMUQQ2c76qRYdWFPWlb3At48YZHclJxYeLKTlSYs76RIuorNG+MUM7ORe7tlNYUWh0KEK0MverRpyqspKV1WlMtu7HYqL/ObPXJ4CTBe1S3yhoy+UnYSYmersRXaleXMBGxpM3LcjoUE6xqHgR8ZHxjI0da3QoncoaU2v6gnZcZBzW/lYpaAtTkUThQ77+oAI7QUz5nnk+FWutWVyymFkps1DK3CvJZU1vaTn+YZnBkXQuNylXhsgKU5FE4UNWrLKgsJM73SQN94DtB7dTXl1u6vqEU8alLQXtz81d0J6cPJmy6jJKj5UaHYoQgCQK39HYyMq9yaQOrKRfP6ODOckn6hMtTha0zXXprj2ZeCfMRhKFj7Bv2MhKew55E82xkp3TouJFJPdNZkT/EUaH0rXgYLIG7jZ9QTs9IZ3QoFApaAvTkEThI7Z9sJUj9GfKheY5nbBrO1+WfMnZKWebvj7h5AsztMOCw8hIyJCCtjANSRQ+YsXntQDkXTzQ4EhOKtpfRFVtlU/UJ5xaW477QEG7oLyAJnuT0aEIIYnCV6ws6suA0GpGjTbPJ3dnfcKXEoXPzNBOzuVE0wk27ttodChCeCdRKKUuUEptU0rtVErd7+LxMUqplUqpeqXUz7uzb0DYv58V1WnkWfdjpis8i4oXYe1vZVjMMKNDcVtU5iifaTkOsjSqMAePv1qUUkHAc8CFwDjgKqVU+wYUh4CfAU/3YF+/d/iLtWxhHFOmm2e0TrO9mSW7l3D2cPOPdjpFcDBZsbspKE0wOpJOpcSkEBsRK4lCmII3PlblADu11jatdQPwFnBp2w201vu11muAxu7uGwhWve9YlSjvskSDIzlpfeV6jtQdYVaK71x2cmqdoV3WbHQoHVJKkZuUKyOfhCl4I1EkAXvb/Fzacl+v7quUmquUKlBKFVRVVfUoULNaudpCEE1MmhZmdCitfLE+4ZQ1rWWGtolbjoPj8tPWA1s5UnfE6FBEgPNGonB1VV339r5a6/la62ytdXZcXJzbwZleUxMr9g5hQlwFUVFGB3PSouJFjIkdQ2K0ec5y3JXxveG+UdBumXi3pmyNwZGIQOeNRFEKDGnzczLg7ke5M9nXLzRv2MRqezZT0s0z0a6xuZFle5b5Xn2iha8UtHOSclAoufwkDOeNV8oaYKRSKkUpFQpcCXzohX39wqb3dlBDNHkXxhgdSquC8gJqGmp8sj4BQFCQTxS0+4X3Y0zsGCloC8N5PFForZuA24BPgS3AO1rrIqXULUqpWwCUUglKqVLgbuBXSqlSpVTfjvb1dMxmsvJzx0pnUy41z+U0Z31i5vCZxgZyBnyhoA2O+RSry1ajtbtXa4XofV4599Zaf6y1HqW1HqG1/k3LfS9qrV9s+b5Sa52ste6rtY5p+f5YR/sGkhWbY4gPO8zwFPNMoFhUvIgJ8ROIjYg1OpQe85UZ2pOTJnOg9gC2wzajQxEBzNwXaQPdwYOsqE5jylnmmWhX31TPV3u/8tn6hFPrDG2TtxzPTW7pJCuXn4SBJFGY2P7/rWUXZ5E3PdjoUFqtKl1FXVOd79YnWpwsaJtnEqMraYPSiAiJkJbjwlCSKExs5Qf7AZhy+WCDIzlpccliLMpC/rB8o0M5Mz5S0A62BJM9OJtVZTLySRhHEoWJrfw6iBDV2DpBzAwWFS8iMzGTmPAYo0M5Yz5T0E7KZX3leuqb6o0ORQQoSRRm1dzMitIhZA4qJTzc6GAcahtrWVW6yufrE05Z+ZGA+QvauUm5NDQ3sL5yvdGhiAAlicKkGjdsZk1zJnnpdUaH0uqrPV/RaG/0+fqEU8alQ31ihrazk6xMvBNGkURhUuv/tZM6+jDlohijQ2m1qHgRwZZgpg2dZnQovcJR0N5u+oJ2Ut8kkqKTZOSTMIwkCpNaucjRsiPvu+Ypti4uWUxOUg5RoSZqOnUmnAXtMvP8jjuSmyydZIVxJFGY1IrNMQwJryJ5iDkmUByrP0ZBeYHf1Cecssb6RkF7ctJkio8UU3XcvzojC98gicKMDh9mZXUqU0aa501h2e5lNOtmv6lPOLXO0P6g1OBIOicT74SRJFGYUNlH69nDMFNNtFtUvIiwoDDykvOMDqVXnZyhfdjoUDqVlZhFkAqSy0/CEJIoTKh1ot0Pkg2O5KTFJYvJG5JHnxDzzOnoDVEZI32ioB0ZGsn4+PFyRiEMIYnChFasCSZc1TExL8LoUAA4dOIQ6yvX+119AvCtgnZSLl+XfY1d240ORQQYSRRmY7ezsnQIk+L3EhpqdDAOS0qWoNF+V59wai1olzYZHUqncpNyOVZ/jK0HthodiggwkihMpm7DNgqb08nLNM9Eu0XFi4gIiSAnKcfoUDwiK9/ZctzcBW3nxDtpECi8TRKFyax9ZyeNhDLlogFGh9Jqcclipg2dRmiQSU5xelnGpb5R0B4dO5p+Yf2kTiG8ThKFyaxwTrQzScfYfTX7KKoq8s/6RIuTBW3zjDJzxaIsTEqaJCOfhNdJojCZlVv6MyKinEHx5pho92XJlwB+W58AHAXtOPO3HAfHxLuN+zdyvOG40aGIACKJwkT0kaOsqE4jb+QBo0Nptah4EX3D+pKZmGl0KB6VNaaW8kYfKGgn52LXdgorCo0ORQQQSRQmsvs/31BJIlNmhBgdSqvFJYvJH5ZPsMXcl2XO1MmW4+YuaOcmOWZoy+Un4U2SKExkxQeOlh15s4cYHIlD6bFSdhza4df1CafWGdqfmbugHRcZh7W/VQrawqskUZjIyoJgoizHSZtsju6si4sXA35en2jRWtDeYO4Z2uA4q5AhssKbJFGYhd3OitKh5CTsIdgkV3kWlSxiQJ8BTIifYHQonmexkBW3h4LSRKMj6VJuUi5l1WWUHjP3ZTLhPyRRmMTx9TvY0JzGlExzrIustWZR8SJmDp+JRQXGn0nWWN8oaMvEO+FtgfEO4APWvFNMM8HkmWSiXfGRYvYc3RMQ9QknX2k5np6QTmhQqNQphNdIojCJFYscLTsmf98cHWMDqT7hlPG94T4xQzssOIz0hHQZ+SS8RhKFSazc2p8xkXsYEGuO/5JFJYuIj4xnbOxYo0Pxmqj0sxhjMX/LcXBMvCsoL6DJbu7LZMI/mONdKcDpY9WOFe1GHTQ6FMBRn1hcvJhZKbNQyhwzxL3CYiErdg8FZT5Q0E7O5UTTCTbu22h0KCIASKIwgR3vbeQgseTlm2Oi3baD26ioqQio+oSTrxS0nRPvpE4hvEEShQms+NBxJjHlyqEGR+IQiPUJp9YZ2u/vNTiSzln7W4mNiJVEIbxCEoUJrCwIJsZylDE5fY0OBXDUJ4b0HcKI/iOMDsWrEhIg/9FzAbj49hSUAqUc95uNUorcpFwpaAuvkERhNK1ZUTaMyYl7sJjgf8Ou7XxZ8mXg1SeAffu6d7/RJidPZuuBrRypO2J0KMLPmeCtKbAdXWejqHkMeVnGTrTLeCkD9Ygi6NdBHKg9wN83/B31iCLjpQxD4xIdc9Yp1pStMTgS4e8kURhs9Zs2NBamXDzQ0DjykvNOW8EuNCiUKclTDIpIdGVS0iRAOskKz5NEYbCVX9ajsJMze5ihcczLn3daq44gFcS8GfMMikh0JSY8hrGxY6WgLTzOJO3nAteKbQMYH1VM3xhjC8eJ0YnMHD6T/+38H+A4m7g+/XoSokxYyRWA43LhlgNb2HJgC+qRk/Wk9IR01t28zsDIhL+RMwoD2auPs6o6lbxRh4wOheMNx9m4byMKxxtOIJ5NxPc70a37jZaXnEeQOnUWuVwuFJ4gicJAm/+1mWP0Y8rM0K439rDHlj5GWXUZl46+FIuyBOTZROWRPuh/vI4eNhytLPw19FYAXv5TncGRuTYvf95pKw8GYoIXnieJwkAr/9Oyot2VxtYnNldt5umVT3Nd+nU8f9HzTBs6LXDfbObMgZISsNu5Zu1dpGDjkXtr0NrowE6XGJ3IdenXtf4slwuFp3glUSilLlBKbVNK7VRK3e/icaWU+mPL498opTLbPFailNqolFqvlCrwRrzesqIgjNigQ5yVHWNYDFprfvrxT4kOjea33/otidGJLLluibzZACGpo/hl3mIK9g3h49fN2VH2oRkPEWJxtH7RWgdughce5fFEoZQKAp4DLgTGAVcppca12+xCYGTLbS7wQrvHZ2mt07XW2Z6O12u0ZmX5UPISd2PkvLbXN77OlyVf8sS3niAuMs64QEzq2r/OYDjFPPLzatOeVdyYcSMAIZYQokLNsYyu8C/eOKPIAXZqrW1a6wbgLeDSdttcCvxdO6wCYpRS5m/heQYOrt3NtuaRTMluMCyGI3VHuGfhPeQk5XBT5k2GxWFmIWPP4pdTvmTNvqH87w3jBx248uCMB0mPT6e2qZYnlj9hdDjCD3kjUSQBbTuslbbc5+42GliolCpUSs3t6CBKqblKqQKlVEFVVVUvhO1Zq94sBiDPwIl2v1r0Kw7UHuDFi14MmOVOe+Lal/MZRgmP/PyYac8q1t2yjqvHX83TK55m95HdRock/Iw33h1cXVhp/3LrbJupWutMHJenfqqUynd1EK31fK11ttY6Oy7O/JdQVnzZQBBNTPpBiiHHLygv4Pk1z3PbpNvISJQ2HZ0JHTuC/5uyhNWVw/n0DXOsGeLKE+c8gUVZ+MXnvzA6FOFnvJEoSoEhbX5OBsrd3UZr7fy6H3gPx6Usn7di20DSo3cREe391dSa7c3c+tGtxEfF8+tZv/b68X3Rda/kM5Tdpj2rABjSbwj3TrmXt4ve5qs9XxkdjvAj3kgUa4CRSqkUpVQocCXwYbttPgSubRn9NBk4qrWuUEpFKqWiAZRSkcB5wCYvxOxRTcdq+bpmLFNGG/PpdH7hfArKC3jmvGfoF97PkBh8TejoFP5v6lJWVabw2ZsHjA6nQ/dNvY+k6CTu+N8d2LXd6HCEn/B4otBaNwG3AZ8CW4B3tNZFSqlblFK3tGz2MWADdgJ/AX7Scn88sFwptQH4GvhIa/0/T8fsad+8s5VaIsmbGe71Y++r2ccDXzzAOSnncGXalV4/vi+7/pXpDGEPD99j3rOKyNBInvjWExRWFPLahteMDkf4CaXN+hd/BrKzs3VBgXmnXDz33c+47YNzKSk8yLBM7xazr33vWt7a9BYbb93I6NjRXj22P3hh2uv85Ks5LPzHfs6dM8jocFyyazt5f81j79G9bL99uwyZFW5RShV2NAVBhroYYEVhKIODKhma4d0ksaRkCa998xr3Tb1PkkQP3fDKNJLZa+pahUVZePb8Z6moqeDJ5U8aHY7wA5IovE1rVlYMJy9pj1cn2jU0N/CTj3/C8Jjh/N/0//Pegf1M2KhhPDBtOV9VnsWiN0269B2QNySPq9Ku4umVMlxWnDlJFF5WWVBKcfMwpmQ3evW4z656ls1Vm/nThX8iIiTCq8f2Nze+MpUkSk1dqwB48ltPolAyXFacMUkUXrbyrRIA8r4T67Vj7j6ym0eWPMJ3x3yXi0dd7LXj+quwkUN5YPpylleOZPFb5j2rkOGyordIovCyFV82EEo9mbO9t1DRnZ/eCcCz5z/rtWP6uxtfmcZgynjk7qOmPqu4b+p9DI4ezJ2f3inDZUWPSaLwspXbB5LVdwdhkd5ZXPC/2//L+1vf58H8BxkWY2w7c38SflYy909fwdLKUXz5ZoXR4XQoMjSSJ855goLyAhkuK3pMEoUXNRw9QUHNGKaM8U5zudrGWm7/5HbGxY3jrry7vHLMQPLjBVNIpJxH7jlmdCidmjNhDjlJOTzwxQPUNNQYHY7wQZIovGjd29upJ9xrE+0eX/Y4JUdKeP7bzxMaZPwqev4mfEQS9+evYEnlaJa82b4rjXnIcFlxpiRReNGKjxxnEnlXe74R4LYD2/jtV7/lmgnXMGP4DI8fL1D9+JUpJFDBI3cfNTqUTslwWXEmJFF40crCMIYFlzJ4ome72zpXrYsMjeSpc5/y6LECXZ8Rg7k/fyWLK8ey9I1So8Pp1BPfegKF4v4vTltkUohOSaLwFq1ZUZHClKQ9Hj/UW5ve4oviL3j87MeJj4r3+PEC3dwFeY6zCpPXKob2G8q9U+7lrU1vyXBZ0S2SKLxk79cVlNkTyZvU5NHjHK07yt0L7yZ7cDZzszpc50n0oj7WRO6bsZpFleNY9sberncwkAyXFT0hicJLVrzpuC485RLPTrR7cPGD7KvZxwsXvUCQxftrXQSqmxfkEc8+059VtB0u+49v/mF0OMJHSKLwkpVLG+lDLRMuH+mxY6yrWMef1/yZW7NvJXuwyyaQwkMiUuK5b8ZqvqhM5as3PX958UzIcFnRXZIovGTF9lhy+m4jJCLEI89v13Zu/ehWYiNi+c05v/HIMUTnbnl1MoPYb/oRUM7hsuXV5TJcVrhFEoUXnDhSz7rjI8kbe9hjx3h57cusLlvN7877HTHhMR47juhYxPBB3Dvjaz6rHN96qdGsZLis6A5JFF5Q8NZOmghhyizPTLSrOl7F/Z/fz8zhM5kzfo5HjiHcc+urucRRZfqzCnAMlwVkuKzokiQKL1jZMtFu8hzPNAK87/P7qG6o5vlvP4/y5iIX4jSRw+O4d+YaFlZOYNWbxUaH06m2w2VX7F1hdDjCxCRReFBCAigFv/jvdAAGjY9HKcf9vWXZ7mW8uv5Vfp73c8bGje29JxY99pNXc4jlgE+cVfxi6i8cw2X/J8NlRcckUXjQvg6WKujofndkvJSBekS13vJfzQfg4x0f9/xJRa+KHBbLvTPX8L/KdFa/aTM6nE45h8uuKV8jw2VFhyRR+Ji85LzTGvwFW4KZNnSaQREJV06eVRwxOpQuzZkwh0mDJ8lwWdEhSRQ+Zl7+PCzq1P+2EEsI82bMMygi4UrUsIHcM7OQTyoz+fqNnUaH0ymLsvDsBY7hsr/96rdGhyNMSBKFh/zpJ1s88rzxUfHkJuW2/hxiCeH69OtJiOrFwofoFT99dRID1UF+fY/5axVThkzhqrSreGrFU+w5au4Jg8L7JFH0suZmuOvcjfzshd4vLG/av4n8Bfks2b2k9awi2BIsZxMmFT1sAPfMXMtHlVmseWOH0eF0yTlc9hef/8LgSITZSKLoRcer7Vw+bjPPfj6eO4a9T/wg16NI4rvZ0PV4w3F+8dkvyHgpg60HtvLqpa8yN3MuFmWRswmTu+3VbAaoQz5xViHDZUVHvLNwcwCo3F3PdzJLWXtoNH+Y8hY/W/J9ng0+8zz80faP+OnHP2X30d3ckH4Dvz33twyMGMh5I85j84HNcjZhctFD+3PPzM/55eJvUfjGNrKuHm10SJ36YOsHAEx9Zeop96cnpLPu5nVGhCRMQM4oekHRiqNMHn2IzYcSeO/a9/nZ8isg+MxycOmxUi5/53IufvNiIkMjWXrdUv566V8ZGDEQgMToRJZct0TOJnzAHzaeDUD2nNEoReutN+fT9JapQ6cSrE792w0NCmVK8hSDIhJmIIniDH3xWjlTpyvq62Hp419xyd8ud7wL9FCTvYlnVz3L2OfG8vGOj3n87MdZd/M6pg+b3otRC2/af8D1y+xM5tN4yrz8eQQHnZooglSQnLkGOEkUZ2DBvF1ccG0cyaqMVf8sJeuB887o+daUrSHnLznc9eldTB86naKfFPHA9AdOmzchhKckRidyffr1hFhOdjkeGDGQTfs3obU2MDJhJEkUPaA1zJu9lRseG8HMPqv5alUww74/qcfPd7TuKLd9fBu5L+ey7/g+/jn7n3x09UdY+1t7MWphRo019UaHcJp5+fNaF70KtgRT11jHua+dS+b8TN7Y+AaNzY0GRyi8TRJFN9XXaX6YvYXH/jWGG+M+5OOdo+mX3bPFiLTWvL3pbcY8N4YXCl7gtpzb2PLTLXx/3PeluV+AGB1TyV+u/IKGI7VGh9LKeVZhURbmZs5l7917efk7L1PXVMecf8/hrD+dxe9X/p7q+mqjQxXeorX2u1tWVpb2hAP7mvT0wTs0aP146j+0/Xhtj59r58Gd+vzXztc8jM56KUuvKVvTi5EKM3Gcg7q+5UQXadA62VKq/3TJQn1i31Gjw9Vaa11+rFznL8jXFdUVrfc125v1h1s/1NNfma55GB3zRIy+/7P7dfmxcgMjFb0FKNAdvKca/qbuiZsnEsWODcf1yMhSHcYJ/eZFr2nd3Nyj56lvqtePLXlMhz8WrqMfj9Z/WPUH3dTc1MvRCjOJj3edJOLjtbbbtf706W/0tJiNGrROUBX6d+d+omt2HzA67E6t2rtKf/+d72vLIxYd+miovuH9G/Tm/ZuNDkucgc4ShXI87l+ys7N1QUFBrz3fiv8c5JLLgqCpiQ/uWsLUZy53a7+MlzJYX7ne5WOzx83m9+f/nqS+Sb0Wp/BdWsOS+dt49Ff1LDowgTiquHv6Gn76ShbRZ3VzhqYX7Tq0i2dWPsOC9Qs40XSCi0ddzM/zfk7+sHy5fOpjlFKFWutsV49JjaILbz9TytmXRNG/+QArX/zG7SQBrju9Alw08iLemf2OJAnRSimYefNovqiawFev7SIrsZwHln2bYSND+HXufzmyca/RIbo0YsAInrvoOfbctYeHZzzMqtJVzPzbTHJfzuWfRf+k2d5sdIiiF8gZRQe0hifn7uKBl0cwNeRr3v8klNhz0t3YT1NypIR1letYsnsJf/76z6csCNMnuA+2O2wyUU50ac2/9/Kbuw7wwZ4M+nKU2yYs4675Y4nN9cxKib2htrGWv2/4O79b+Tt2HtqJtb+Vuyffzfy18/lm3zenbS8zvs2jszMKSRQuNDbCT87dwctLRnJl349Z8HUq4aOHnbZdk72JbQe2sa5yHWsr1rKuch3rK9dzpO4I4Jio1DesL0frjmLHTmhQKDdl3MRzFz3X49hE4NnwaSW/ua2Cf+2cSAS13DpmMfc8fxbpV411OWkvPh4qK70fZ1vN9mY+2PYBT614ilWlqwgLCqPJ3kSzPnmGIa8Hc5FE0YWQ2zJoilt/2v2qYgJNjy/BMiCGuqY6Nu7b2JoQ1lWu45t931DXVAdAeHA4E+InkJGQ4bglZjB+0HiO1B3B+kcrdU11cjYhzsjmpQd4/NY9vLl5IqE0UEefDrc1y8taa81Xe7/i0aWPsnDXwlMeC7WE8t4V75GTnMPAPgOlpmEwSRRdUBf/BDL+CsENJ+9sCoayXK75jpV1levYUrWl9dNQv7B+ZCRmnJIUxsSOIdjiur/TTz76CS8VvsQtWbfIpydxxnasOcL/m1vMgvUZHW6jDx6C/v273U4mIcF1a5HeOEu56t2reKfoHZdrc/cN64u1v5UR/Ucwov8Ix/cDHN8P6Tekw9dWRwNG5JJW9xmeKJRSFwB/AIKAl7XWT7R7XLU8/m2gFrhOa73WnX1d6XaiiK6AO6wQUnfaY4lRiaclhZSYlG59+qmoruDKd6/k7e+/LWcTotd09icYw2ESVSWJYYcYHHWMxP51JA6yk5hkYfDwUBJHRpE4rj9RZyVAbCxYLF0+55m8VSQkwL7jbV5njX3gH58QE3+Mh/9gY9fhXY7boV0UHymmofnkh7ZgSzDDY4afkkhGDHAkk9yf/Ym6UX9v9yEvlD5bb6L2nz3/UOaJBOSppNZbz9tZovB4m3GlVBDwHHAuUAqsUUp9qLXe3GazC4GRLbdc4AUg1819z1xNIqy/HrJfBKWhOQh2XgD/eZny6jN/Y3d2ehXCW344bTfl+4OoODSI5TVWKg4OoH5H2GnbRVHNYHaQGHqIxKhjwPkdPueGPy4hOCyIkD7BhPQJJrhPiOP7iBDH9xGOW3BEKCosFEJDHV2UlWo5S2l5nWW9BOuuh90zOLIb7ph86nGa7c2UV5e3Jo62SWRN2RoO1x0+ufE4oH0CszRzorGWny/8OX2C+9AnpM8pX8ODw0+7r/3XnKQcNldtPiVhnWkX3bzkvF5/Tk8+b1seP6NQSuUBD2utz2/5+QEArfX/a7PNS8CXWus3W37eBswEhne1ryvdPqNQQFS7Tzt/sEFNgmmu9QrRXnc+/WsNhw9DRUk9FZsPU7G9moriOspLm6nYZ6HiUBgV1VHsrO2dIdtBNBFCI8EtXw8zwPFAVAV8/0r419tQ4/gQNmoUp7Rf7+rWHHqYhshdNETuwnZkF0x4DQZuA4t2JI2GSGgOo0/fE5xoOtEr/x6n/uH9CbYEY1GWU25KqdPua30Mx2PN9ma2HNiCbpPZFIqJ8RMJCTrZhNHV1QqF6/9spRQNzQ2sq1h3yvP2pB5q6BkFkAS0HQReiuOsoattktzcFwCl1FxgLsDQoUO7H2VNu087NXKJSJhbfHzH9YT2lIIBA2DAgDBSMxMA13/fnSWffz+1i8a6Zprqm2msa6ax3k5TQzON9brlezuNDdpxa9Q0NWoaGxyjCP+8tuXTbU0ivHrq2XVGRmdNTlzd+qN1NropG9syHK/XO6xgqYOmPvCnnVCTQG1L54mG5gZONJ3gROMJl1/rmupcPvbvrf9mQ+UGmnUzFmVhXOw48oflo9HYtd3lraPHtHbcf7zxOHuO7kGjUSiG9B1CYnRi6+9Cn3Z65Pg3uNJ226S+SZQdK0OjCQ0K7fWVL72RKFz96bX/l3e0jTv7Ou7Uej4wHxxnFN0JsPUFt2QexBXB0nmt9wthVt4eAvu9n/d8/safO0lAb73V46d1JLZOPuQppQgLDiMsOIyY8JhuPfcNGTdg/aOV5qZmwoLC+Ozaz874zbeiuqJ1FGR4cDirf7y6V97Q2z6vJ9YP8cbM7FJgSJufk4FyN7dxZ98zVlnZ8kmlOhG9YAm62nHJyeix6EJ4W0cfjkz/oWnJPNgzrfVDXm9o20W3tz6he+I5Pfm8rTpqAtVbNxxnLTYgBQgFNgCp7ba5CPgExxnEZOBrd/d1dfNU91ghRM901hjRjM/r5KqLrhmfszeeF6ObAiqlvg08i2OI6yta698opW5pSVQvtgyP/TNwAY7hsddrrQs62rer4/V2U0AhhPB3hs+j8DZJFEII0T3SPVYIIUSPSaIQQgjRKUkUQgghOiWJQgghRKckUQghhOiUJAohhBCdkkQhhBCiU5IohBBCdEoShRBCiE755cxspVQVsNvoONqJBQ4YHYSbJFbP8aV4fSlW8K14zRjrMK11nKsH/DJRmJFSqqCj6fFmI7F6ji/F60uxgm/F60uxglx6EkII0QVJFEIIITolicJ75hsdQDdIrJ7jS/H6UqzgW/H6UqxSoxBCCNE5OaMQQgjRKUkUHqSUGqKUWqyU2qKUKlJK3WF0TF1RSgUppdYppf5rdCxdUUrFKKX+pZTa2vI7zjM6po4ope5q+RvYpJR6UykVbnRMbSmlXlFK7VdKbWpz3wCl1GdKqR0tX/sbGWNbHcT7VMvfwjdKqfeUUjEGhtjKVaxtHvu5UkorpWKNiM1dkig8qwm4R2s9Fsda4D9VSo0zOKau3AFsMToIN/0B+J/WegwwEZPGrZRKAn4GZGut03As63ulsVGd5lUcSxG3dT/whdZ6JPBFy89m8Sqnx/sZkKa1ngBsBx7wdlAdeJXTY0UpNQQ4F9jj7YC6SxKFB2mtK7TWa1u+r8bxRpZkbFQdU0olAxcBLxsdS1eUUn2BfOCvAFrrBq31EUOD6lww0EcpFQxEAOUGx3MKrfVS4FC7uy8F/tby/d+A73ozps64ildrvVBr3dTy4yog2euBudDB7xbg98B9gOkLxZIovEQpNRzIAFYbHEpnnsXxh2s3OA53WIEqYEHLpbKXlVKRRgflita6DHgaxyfHCuCo1nqhsVG5JV5rXQGODz3AIIPj6Y4bgE+MDqIjSqlLgDKt9QajY3GHJAovUEpFAe8Cd2qtjxkdjytKqYuB/VrrQqNjcVMwkAm8oLXOAI5jrksjrVqu7V8KpACDgUil1A+Njcp/KaV+ieOy7+tGx+KKUioC+CXwoNGxuEsShYcppUJwJInXtdb/NjqeTkwFLlFKlQBvAWcrpf5hbEidKgVKtdbOM7R/4UgcZvQtoFhrXaW1bgT+DUwxOCZ37FNKJQK0fN1vcDxdUkr9CLgYmKPNO/Z/BI4PDRtaXm/JwFqlVIKhUXVCEoUHKaUUjmvoW7TWzxgdT2e01g9orZO11sNxFFoXaa1N+6lXa10J7FVKjW656xxgs4EhdWYPMFkpFdHyN3EOJi28t/Mh8KOW738EfGBgLF1SSl0A/AK4RGtda3Q8HdFab9RaD9JaD295vZUCmS1/06YkicKzpgLX4Ph0vr7l9m2jg/IjtwOvK6W+AdKBx40Nx7WWs55/AWuBjThed6aamauUehNYCYxWSpUqpW4EngDOVUrtwDE65wkjY2yrg3j/DEQDn7W81l40NMgWHcTqU2RmthBCiE7JGYUQQohOSaIQQgjRKUkUQgghOiWJQgghRKckUQghhOiUJAohhBCdkkQhhBCiU5IohDCAUuqvSqmLjI5DCHdIohDCGOmAT3QOFSLY6ACECARKqVHAK0C/lq8JWutSY6MSwj1yRiGEh7UsVvQP4G6t9XhgJLDV2KiEcJ+cUQjheZfh6CD8dcvPRcAJA+MRolvkjEIIz5sAtF0QKgupTwgfIolCCM87CKQBKKWygKuA9UYGJER3SJtxITxMKRULfAyEAtuAfGBoy2p3QpieJAohhBCdkktPQgghOiWJQgghRKckUQghhOiUJAohhBCdkkQhhBCiU5IohBBCdEoShRBCiE5JohBCCNGp/w9od1q1sCxv9AAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(6,6))\n",
    "plt.plot(all_sp1, Pl, 'ro-', label='original network')\n",
    "plt.plot(all_sp1_d, Pl_d, 'bs-', label='degree-preserving random network')\n",
    "plt.plot(all_sp1_r, Pl_r, 'gv-', label='random network')\n",
    "plt.legend(loc=0)\n",
    "plt.xlabel(\"$d$\")\n",
    "plt.ylabel(\"$p_d$\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "d405888a-533c-42bb-acf0-adcd9011bdb0",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "True\n",
      "False\n"
     ]
    }
   ],
   "source": [
    "print(nx.is_connected(LCC))\n",
    "print(nx.is_connected(G_d))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "a93490e1-7bb8-4fc6-8c19-4683f33919ad",
   "metadata": {},
   "outputs": [],
   "source": [
    "k1, Pk1 = get_pdf(LCC)\n",
    "k2, Pk2 = get_pdf(G_d)\n",
    "k3, Pk3 = get_pdf(G_r)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "c4b36c12-485a-4a2e-902e-1e52e27cfa25",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAE4CAYAAABMlxQUAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAZKElEQVR4nO3df2zc933f8ef76Jwo1pbWOq5JixYt1kFQLYWkRcosz5AGLC0cIG4SJ2gd37qa88JIQrr8s2HNPAKtPc4ZtGGrASGy8sOep5tjw1HTJFsRYAMiw5s9RK41JE7Q2qAcRTYlOW0h2bKYRNFnf1CWKVFH3vd43/veffl8AAeZn+/3vt83mS/5yvc+3+/7GyklJElqVqXoAiRJvcXgkCRlYnBIkjIxOCRJmRgckqRMDA5JUiZXFV1Ant797nenm266qegyJKmnPP/88z9JKV3XaHkpgyMi7gDuuPnmmzl06FDR5UhST4mIHy20vJQfVaWUvplSGl+9enXRpUhS6ZQyOCRJ+TE4JEmZGBySpEwMDklSJgaHJCmTUgZHRNwREftOnTpVdCmSVDqlDI4lXY5brzM9vIXtcZDjN26Ber39BUpSDytlcLSsXofxcR54dYxnuI37j43B+LjhIUlzRJmfALh58+aU5c7xlZUZZlL/vPH+mOHs+fnjklRGEfF8Smlzo+WeccwxlUa5mzoDnAFggDPU2M+RtK7gyiSpexgccwyNVFnFaWbop5+zzNDPKk4zOLKi6NIkqWsYHHNNTnKicgM72Mtz3MIO9nK8sgYmJ4uuTJK6Rim747asVuMAdbhvNxw9yp61u2dDo1YrujJJ6hqlnByf01b9Uy+99FLR5UhST1mWk+O2VZek/JQyOCRJ+TE4JEmZGBySpEwMDklSJgaHJCkTg0OSlInBIUnKxOCQJGVSyuDwCYCSlJ9SBod3jktSfkoZHJKk/BgckqRMDA5JUiYGhyQpE4NDkpSJwSFJysTgkCRl4jPHS2LTw5s4fPzwvPGNgxt54dMvdL4gSaVlcJTE1uGt/OD49/kZ5y6OVbmKW4dvzWV/BpW0fPlRVUlM/GQ9lZ+fu2Ss7+fnmPjJ+lz2t3V4K9W+6iVj1b5qbkElqXsYHCUxOvbPmDm8E85d+GN+rsrZF3axbuzeXPY3sW2CSlx6+PRFHxPbJ3LZn6TuYXCUxFQa5WMH10O68D9pqnDn07/OkbQul/0NXTPE2Maxi2cd1b4qYxvHGLx6MJf9SeoeBkdJDI1Uuf7NPjj8+3C+Ai/cw/VvVhgcWZHbPueedXi2IS0fpQyOZdlWfXKSE5Ub+P2Da/l7R1dyz9NrOV5ZA5OTue3y7bOOSlQ825CWkUgpFV1DbjZv3pwOHTpUdBmdU6/DfffB0aOwdu1saNRque5y+o1p7vraXTzxiScMDqkkIuL5lNLmhssNDknSXIsFh/dxqHDeEyL1llLOcai3eE+I1FsMDhXOe0Kk3mJwqHDeEyL1Fuc45vCz9uJMbJvgkcOPAJ5tSN3OM445/Ky9ON4TIvUOg2MOP2sv1sS2CW5be5s/b6nLGRxz+Fl7sYauGeLgPQf9eUtdzuC4jP2XJGlhBsdl/KxdkhbmVVVXMLFtghdff9GzDUm6AoPjCt7+rF2SNJ8fVUmSMjE4JEmZGBySpEwMDrWuXmd6eAvb4yDHb9wy+yApSaVncKg19TqMj/PAq2M8w23cf2wMxscND2kZ8AmAasnKygwzqX/eeH/McPb8/HFJvWOxJwD2zBlHRIxGxJcj4qmiaxFMpVHups4AZwAY4Aw19nMkrSu4Mkl568h9HBHxFeDDwMmU0vvmjN8O/AnQB3wppfT5RttIKU0B9xoc3WFopMqqH51mhn76OcsM/aziNIMjKzpah63wpc7r1BnHo8Dtcwciog/YA3wIWA98MiLWR8RvRMS3Lnv9aofqVLMmJzlRuYEd7OU5bmEHezleWQOTkx0tw1b4Uud1JDhSSk8Df3PZ8AeAl1NKUymlnwFfBT6SUvpeSunDl71OdqJOZVCrceCxN9kzspsN8T32jOzmwGNvQq3W0TJshS91XpFzHGuAH8/5+tiFsSuKiGsjYi+wKSI+t8B64xFxKCIOvf766+2rVvPVavDKK3D+/Oy/HQ4NsBW+VIQigyOuMNbwEq+U0l+nlHaklH4tpfTgAuvtSyltTiltvu6669pSqLqbrfClzioyOI4BN875ehh4raBa1MNshS91VpHdcb8LvCci1gGvAncBdxdYj3qYrfClzunIGUdEPA48C7w3Io5FxL0ppXPAZ4BvAz8Enkwpvdim/d0REftOnTrVjs2pB/jYWalzvHNcknSJ0tw5LknqDgaHJCmTUgaHcxw9yBbtUs8oZXCklL6ZUhpfvXp10aWoGbZol3qKk+MqnC3ape7i5Li6ni3apd5icKhwQyNVVlF8i3ZJzSnyzvHcRMQdwB0333xz0aWoGZOTnPgnV7Pj/F7G2cc+xpmuDOfWot1neEhLU8ozDifHe0yHW7T7DA9paZwc17Iz/cY0ow+NMnNu5uLYyqtWMvXZKVuWSDg5Ls3jMzykpTE4tCz5DA+pdaUMDu8c12J8hofUulIGh5PjasbEtgluW3ubZxtSRqW8HFdqxtvP8JCUTSnPOCRJ+TE4JEmZGBySpEwMDklSJqUMDi/HlaT8lDI4vBxXkvJTyuCQJOXH4JAkZWJwSJIyMTgkSZkYHOp99TrTw1vYHgc5fuMWqNeLrkgqNYNDva1eh/FxHnh1jGe4jfuPjcH4uOEh5aiUweF9HMvHyt/7OPHWGb7ALs7TxxfYRbx1hpW/9/GF3+hZitSyUgaH93EsH1NplLupM8AZAAY4Q439HEnrGr/JsxRpSUoZHFo+hkaqrOI0M/TTz1lm6GcVpxkcWdHwPS2fpUgCDA71uslJTlRuYAd7eY5b2MFejlfWwORkw7e0dJYi6SKDQ72tVuPAY2+yZ2Q3G+J77BnZzYHH3oRareFbWjlLkfQOg0O9r1aDV16B8+dn/10gNICWzlIkvcNHx2r5qdU4QB3u2w1Hj7Jn7e7Z0FgscCQBBoeWq1pt0aDY9PAmDh8/PG984+BGXvj0CzkVJnU/P6qSGtg6vJVqX/WSsWpflVuHby2oIqk7GBxSAxPbJqjEpb8ifdHHxPaJgiqSukMpg8M7x9UOQ9cMMbZx7OJZR7WvytjGMQavHiy4MqlYpQwO7xxXu8w96/BsQ5pVyuCQ2uXts45KVDzbkC7wqippERPbJnjx9Rc925AuMDikRQxdM8TBew4WXYbUNfyoSpKUicEhScrE4JAkZWJwSJIyMTgkSZkYHJKkTAwOSVImBockKRODQ5KUicEhScqklMFhW3VJyk8pg8O26mqbep3p4S1sj4Mcv3EL1OtFVyQVLlNwRER18bWkkqjXYXycB14d4xlu4/5jYzA+bnho2YuUUvMrRzyYUvrchf/+Byml/51bZW2wefPmdOjQoaLLUI9aWZlhJvXPG++PGc6enz8ulUVEPJ9S2txoedaPqr4dEeMRUQM+tLTSpO42lUa5mzoDnAFggDPU2M+RtK7gyqRiNf08joj4MnAK2Ag8l1L613kVJXWDoZEqq350mhn66ecsM/SzitMMjqwoujSpUE2fcaSU7gXuAyaAkxHxcG5VSd1gcpITlRvYwV6e4xZ2sJfjlTUwOVl0ZVKhFj3jiIgPAX8E/B3g/wH/KaX0n3OtSuoGtRoHqMN9u+HoUfas3T0bGrVa0ZVJhVp0cjwijgD/GPgB8H5mQ2RPSunx3KtbIifHJSm7xSbHm5njODHn6qn/GRHPAv8X6PrgkCS1XzNzHK9ExL+dcw/Hz4E3cqxJktTFmgmOBNwJ/DgingFeBr4TEe/JtTJJUlda9KOqlNInASKiH3gfsOHC60sRMZpSujHfEiVJ3STL5bgzKaVDKaUvp5T+eUppu6EhYT8rLTulbHIodYz9rLQMZepV1Wu8HFd5s5+VyqjdvaokzWE/Ky1HBoe0BEMjVVZhPystLwaHtBT2s9Iy1HR3XElXYD8rLUM9c8YRER+NiC9GxJ9FxG8VXY90Ua0Gr7wC58/P/nt5aHi5rkqmI8EREV+JiJMR8f3Lxm+PiL+MiJcj4g8X2kZK6esppU8B9wC/m2O5Uvt4ua5KqCOX40bENuBN4LGU0vsujPUBfwX8JnAM+C7wSaAPePCyTfzTlNLJC+/7j0A9pfQXi+3Xy3FVNC/XVS9qR3fcJUspPR0RN102/AHg5ZTSFEBEfBX4SErpQeDDl28jIgL4PPDnzYSG1AmbHt7E4eOH541vHNzIC59+gak0yr9gN1/no7zFLzHAGT7Gn/If0r8Epjter9QORc5xrAF+POfrYxfGGvkD4IPAJyJiR6OVLjwT/VBEHHr99dfbU6nUwNbhrVT7qpeMVfuq3Dp8K+DluiqnIoMjrjDW8HOzlNJDKaX3p5R2pJT2LrDevpTS5pTS5uuuu64thUqNTGyboBKX/hr1RR8T2ydmv/ByXZVQkcFxDJjbJHEYeK2gWqSWDF0zxNjGsYtnHdW+KmMbxxi8enB2hVqNA4+9yZ6R3WyI77FnZDcHHnvTy3XV0zrWq+rCHMe35kyOX8Xs5Pg/Al5ldnL87pTSi+3ap5Pj6oTpN6YZfWiUmXMzrLxqJVOfnXonOKQe1BW9qiLiceBZ4L0RcSwi7k0pnQM+A3wb+CHwZLtCIyLuiIh9p06dasfmpAW9fdZRicqlZxtSSdkdV2qD6Temuetrd/HEJ54wONTzuuJyXKnshq4Z4uA9B4suQ+qInmk5IknqDqUMDuc4JCk/pQyOlNI3U0rjq1evLroUaWE2QFQPKmVwSD3BBojqUV5VJRXEBojqVl1xH4ek+XxeuXpVKYPDyXH1AhsgqleVMjicHFdPsAGiepQ3AEpF8Xnl6lEGh1SkWs2gUM8p5UdVkqT8GBySpExKGRxeVSVJ+SllcHhVlSTlp5TBIUnKj8EhScrE4JC6VaPOuXbUVcEMDqkbNeqcu2uXHXVVuFJ2x42IO4A7br755k+99NJLRZcjZdaocy4kIOaN2lFX7bQsu+N6VZV6XaPOuYfZYEddFa6UwSH1ukadczf0/cCOuiqcwSF1o0adc8fH7airwpVyjuNtPgFQPa1eh/vug6NHYe3adzrnNhqX2mSxOQ6DQ5J0iWU5OS5Jyo/BIUnKpJTBYXdcScpPKYPD+zi07NiGRB1UyuCQlpVG7UkMD+XEq6qkHteoPYltSNQqr6qSSq5RexLbkCgvBofU4xq1J7ENifJicEi9rlF7EtuQKCdXFV2ApCWq1ThAHe7bDUePsmftbtuQKFcGh1QGtZpBoY7xoypJUiYGhyQpk1IGhy1HJCk/pQwOW45IUn5KGRySpPwYHJKkTAwOSVImBockKRODQ5KUicEhScrE4JAkZWJwSJIyMTgkSZkYHJKkTAwOSVImBockKZNSBofdcSUpP6UMDrvjSkC9zvTwFrbHQY7fuAXq9ezrLLS8me2rnFJKpX29//3vT9KytH9/SgMDaSd7UoVzaSd7UhoYmB1vdp2FljezffUs4FBa4G9rzK5TTps3b06HDh0qugyp41ZWZphJ/fPG+2OGs+f7m1pnoeXAottX74qI51NKmxstL+VHVdJyN5VGuZs6A5wBYIAz1NjPkbSu6XUWWt7M9lVeBodUQkMjVVZxmhn66ecsM/SzitMMjqxoep2FljezfZWXwSGV0eQkJyo3sIO9PMct7GAvxytrYHKy+XUWWt7M9lVeC02A9PrLyXEta/v3pzQyklLE7L9XmrhebJ2FljezffUknBx3clySsnByXJLUVgaHJCkTg0OSlInBIUnKxOCQJGVicEiSMjE4pBxtengT8ccx77Xp4U1Fl9aapXbEzdqNd9cuO/B2IYNDytHW4a1U+6qXjFX7qtw6fGtBFS1BvQ7j4zzw6hjPcBv3HxuD8fHm/5gv9v4rLf/CF1rfn3LjDYBSjqbfmGb0oVFmzs1cHFt51UqmPjvF4NWDBVaWXTMdd5fy/kbLW92fWucNgFKBhq4ZYmzj2MWzjmpflbGNYz0XGtBcx92lvP/y5St5i5s4wkreaml/yo/BIeVsYtsElZj9VeuLPia2TxRcUWuW2hE3azfen7KCAc7wU1bYgbfLGBxSzt4+66hEpWfPNoCld8RtoRvv3/LLduDtRgt1QOymF/DrwF7gKWBnM++xO666xWunX0vbHtmWpt+YLrqUpVlqR9ys3Xh37rQDbwHohu64EfEV4MPAyZTS++aM3w78CdAHfCml9PkmtlUBvphSunexdZ0cl6TsumVy/FHg9rkDEdEH7AE+BKwHPhkR6yPiNyLiW5e9fvXCe34beAb4Xx2qW5J0mas6sZOU0tMRcdNlwx8AXk4pTQFExFeBj6SUHmT27ORK2/kG8I2I+O/Af8uxZElSAx0JjgbWAD+e8/Ux4O83Wjki/iFwJ7AC+B8LrDcOjAOsXbu2DWVKkuYqMjjiCmMNJ1xSSt8BvrPYRlNK+4B9MDvH0WJtkqQGirwc9xhw45yvh4HXCqpFktSkIoPju8B7ImJdRFSBu4BvFFiPJKkJHQmOiHgceBZ4b0Qci4h7U0rngM8A3wZ+CDyZUnqxTfu7IyL2nTp1qh2bkyTNYZNDSdIluuU+DklSSRgckqRMShkcznFIUn5KGRwppW+mlMZXr15ddCmSVDqlDA5JUn4MDklSJgaHJCmTUgaHk+OSlJ9SBoeT45KUn1IGh9QLNj28ifjjmPfa9PCm/HZarzM9vIXtcZDjN26Bej3b8la228o+mx3L43vetSufn1G71OtM/8rf5ZZ4lq3xLMd/Zf07NXSqtoWeK9vrL585rm6281s7U/WBauKPuPiqPlBNu761K58d7t+f0sBA2smeVOFc2smelAYG3nmO92LLW9luK/t817tSqlYXH1tqbY2WQ/t/Ru2yf39K73pX2smeBL9IcH62hmp19vnsbaqNbnjmeFHsVaVuNv3GNKMPjTJzbubi2MqrVjL12SkGrx5s+/5WVmaYSf3zxvtjhrPn+xdd3sp2gZb22ayl1JZl/0v9GbXLwvUmrvSYo1Zqs1eV1KWGrhlibOMY1b4qANW+KmMbx3IJDYCpNMrd1BngDAADnKHGfo6kdU0tb2W7rezzTr7Gx/jaomNLre1Ky1fyFjdxhJW81dafUbtMpVE+xtfo4+cXx/o4x508xWE2dKy2UgaHV1WpV0xsm6ASs7+GfdHHxPaJ3PY1NFJlFaeZoZ9+zjJDP6s4zeDIiqaWt7LdVvZ5PSe4npOLji21tist/ykrGOAMP2VFW39G7TI0UuV6TvIL+pg9w0j8gj6u5yQb+n7QsdpKGRzJq6rUI94+66hEJdezDQAmJzlRuYEd7OU5bmEHezleWQOTk80tb2W7rewzhjgRQ4uPLbW2Bsv/ll9u/8+oXSYnORGDrOMIv8MT/A5PsI4pjscNMD7eudoWmgDp9ZeT4+oFr51+LW17ZFuafmM6/53t35/SyEhKEbP/Xj5xutjyVrbbyj6bHcvje965M5+fUbvs35/StdemBLOva6+9dPK+DbXh5LiT45KUhZPjkqS2MjgkSZmUMji8qkqS8lPK4EheVSVJuSllcEiS8mNwSJIyMTgkSZmU+j6OiDgFvLTAKquBRjPo7wZ+0vaiOmOh76vb97eUbbXy3mbf08x6i63j8dZ9+/N4u7KRlNJ1DZcudHdgr7+Afa0uZ5E7J7v5tdj33c37W8q2Wnlvs+9pZj2Pt97bn8dba6+yf1T1zSUu71Wd/r7aub+lbKuV9zb7nmbW83jrvf15vLWg1B9VLUVEHEoL3HIvtZPHmzppqcdb2c84lmJf0QVoWfF4Uyct6XjzjEOSlIlnHJKkTAwOSVImBockKRODo0kR8UsR8V8i4osRUSu6HpVbRIxGxJcj4qmia1H5RcRHL/xt+7OI+K3F1l/WwRERX4mIkxHx/cvGb4+Iv4yIlyPiDy8M3wk8lVL6FPDbHS9WPS/L8ZZSmkop3VtMpSqDjMfb1y/8bbsH+N3Ftr2sgwN4FLh97kBE9AF7gA8B64FPRsR6YBj48YXVftHBGlUej9L88SYt1aNkP97+zYXlC1rWwZFSehr4m8uGPwC8fOH/8f0M+CrwEeAYs+EBy/znptZkPN6kJclyvMWsfw/8eUrpLxbbtn8A51vDO2cWMBsYa4ADwMcj4guUt3WEOu+Kx1tEXBsRe4FNEfG5YkpTCTX6+/YHwAeBT0TEjsU2clU+tfW0uMJYSimdAcY6XYxKr9Hx9tfAor/AUkaNjreHgIea3YhnHPMdA26c8/Uw8FpBtaj8PN7USW053gyO+b4LvCci1kVEFbgL+EbBNam8PN7USW053pZ1cETE48CzwHsj4lhE3JtSOgd8Bvg28EPgyZTSi0XWqXLweFMn5Xm82eRQkpTJsj7jkCRlZ3BIkjIxOCRJmRgckqRMDA5JUiYGhyQpE4NDkpSJwSFJysTgkDooIj4YEf+16DqkpTA4pM7aALxQdBHSUhgcUmdtAF6IiBUR8WhE/LuIuFKra6lr+TwOqbM2ACeZbTL3pZTS/oLrkTKzyaHUIRHxLuAnwI+AT6eUni24JKklflQldc56Zp+HcA74RcG1SC0zOKTO2QD8H2YfnvNIRFxfcD1SSwwOqXM2AN9PKf0V8K+AJy98fCX1FOc4JEmZeMYhScrE4JAkZWJwSJIyMTgkSZkYHJKkTAwOSVImBockKRODQ5KUyf8HGIpiYAa9/w4AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x345.6 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(6, 4.8))\n",
    "plt.plot(k1, Pk1, 'ro')\n",
    "plt.plot(k2, Pk2, 'b*')\n",
    "plt.plot(k3, Pk3, 'gv')\n",
    "plt.xlabel(\"$k$\")\n",
    "plt.ylabel(\"$p_k$\")\n",
    "plt.xscale(\"log\")\n",
    "plt.yscale(\"log\")"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
